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Abstract 

Ki 4 , decays offer several reasons of interest: they allow an accurate measurement of TTTr-scattering lengths; they 
provide the best source for the determination of some low-energy constants of xPTj one form factor is directly related 
to the chiral anomaly, which can be measured here. We present a dispersive treatment of Kn decays that provides a 
resummation of tttt- and ATr-rescattering effects. The free parameters of the dispersion relation are fitted to the data of 
the high-statistics experiments E865 and NA48/2. The matching to xPT at NLO and NNLO enables us to determine the 
LECs Li, L 2 and L 3 . With recently published data from NA48/2, the LEG Lg can be determined as well. In contrast 
to a pure chiral treatment, the dispersion relation describes the observed curvature of one of the form factors, which we 
understand as a rescattering effect beyond NNLO. 
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1 Introduction 


Ki 4 denotes the semileptonic decay of a kaon into two pions and a lepton pair. Its amplitude has a similar 
structure to that of i^Tr scattering, with the difference that in K (4 decays one of the axial currents couples to an 
external field, the W boson, which decays into the lepton pair - the of this axial current is therefore variable 
rather than being stuck at as in Ktt scattering. This difference has the important consequence that in 
K (4 decays the allowed kinematical region reaches down to lower energies, E < Mk, whereas in Ktt scattering 
E > Mk + From the point of view of chiral perturbation theory (yPT) [1, 2, 3], the low-energy effective 
theory of QCD, Ki 4 decays offer similar information as Kt: scattering, but in a kinematical region where the 
chiral expansion is more reliable. 

Due to its two-pion final state, K 14 is also one of the cleanest sources of information on tttt interaction 
[4, 5, 6 ]. 

The latest high-statistics K 14 experiments E865 at BNL [7, 8 ] and NA48/2 at CERN [ 6 , 9] have achieved 
an impressive accuracy. The statistical errors of the S-wave of one form factor reach in both experiments the 
sub-percent level. Matching this precision requires a theoretical treatment beyond one-loop order in the chiral 
expansion. A first treatment beyond one loop, based on dispersion relations, was already done twenty years ago 
[10]. The full two-loop calculation became available in 2000 [11]. However, as we will show below, even at two 
loops yPT is not able to predict the curvature of one of the form factors. 

Here, we present a new dispersive treatment of K 14 decays. The form of the dispersion relation we solve is 
not exact, but relies on an assumption (absence of D- and higher wave contributions to discontinuities) that 
is violated only starting at 0(p®) in the chiral expansion. It resums two-particle rescattering effects, which 
we expect to be the most important contribution beyond two loops. Indeed, we observe that the dispersive 
description is able to reproduce the curvature of the form factor. 

The dispersion relation is parametrised by subtraction constants, which are not constrained by unitarity. 
These have to be determined by theoretical input or by a fit to data. It turns out that the available data does 
not constrain all the subtraction constants to a sufficient precision. Therefore, we use the soft-pion theorem, a 
low-energy theorem for K 14 that receives only SU{2) chiral corrections, as well as some chiral input to constrain 
the parameters that are not well determined from data alone. 

The present treatment of A '^4 decays represents an extension and a major improvement of our previous 
dispersive framework [12, 13, 14]. The modifications and improvements concern the following aspects: 

• Instead of a single linear combination of form factors, now we describe the two form factors E and G 
simultaneously. This allows us to include more experimental data in the fits. 

• The new framework is valid also for non-vanishing invariant energies of the lepton pair. In the previous 
treatment, we neglected the dependence on this kinematic variable. This approximation is no longer used 
and the observed dependence on the lepton invariant energy can be taken into account. 

• We apply corrections for isospin-breaking effects in the fitted data that have not been taken into account 
in the experimental analysis. 

• We perform the matching to yPT directly on the level of the subtraction constants, which avoids the 
mixing with the treatment of rescattering effects. 

• Besides a matching to one-loop yPT, we also study the matching at two-loop level. 

The first two points required a substantial modification and extension of the dispersive framework from the 
very start, but rendered it much more powerful. The old treatment can be understood as a limiting case of the 
new framework. 

The outline is as follows: in section 2, we derive the dispersion relation for the Ki 4 form factors, which has 
the form of a set of coupled integral equations. In section 3, we describe the numerical procedure that is used 
to solve this system. Section 4 is devoted to the determination of the free parameters of the dispersion relation 
and the derivation of matching equations to yPT. In section 5, we present the results of the fit to data and the 
values of the low-energy constants L\, L 2 and Lg obtained in the matching to yPT. Section 6 concludes the 
main text. The appendices contain several details on the kinematics, the derivation of the dispersion relation 
and explicit expressions for the matching equations. Further details that are omitted here can be found in [15]. 
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2 Dispersion Relation for K£ 4 ^ 

2.1 Decay Amplitude and Form Factors 

Kii are semileptonic decays of a kaon into two pions and a lepton-neutrino pair: 

K+{k) (1) 

where (. € {e,/r} is either an electron or a muon. There exist other decay modes involving neutral mesons. 
Their amplitudes are related to the above decay by isospin symmetry - in our dispersive treatment of we 
will work in the isospin limit and could therefore describe as well the neutral mode. In the present analysis, 
however, we only consider the charged mode because it is the one which has been measured more accurately. 

In the standard model, semileptonic decays are mediated by W bosons. After integrating out the W boson 
from the standard model Lagrangian, we end up with a Fermi type effective current-current interaction. The 
matrix element of Ki 4 ^ then splits up into a leptonic times a hadronic part. The leptonic matrix element can be 
treated in a standard way. The hadronic matrix element exhibits the usual V — A structure of weak interaction: 

out(7r+(pi)7r“(p2)^’^(w)t'£(Piy)|Ar+(fc))in = i{2TT)'^5^‘^\k - Pi - P2 -pi-pu)T, (2) 

T = ^V*,u{p^)-i^^{l - 75)i’(w)(7r+(pi)7r-(p2)|V^(0) - A^{Q)\K+{k)), (3) 

where and Note that although we drop the corresponding labels, the meson states are 

still in- and out-states with respect to the strong interaction. 

The Lorentz structure of the currents allows us to write the two hadronic matrix elements as 

V;- := {7r+{pi)n-{p,)\V^{0)\K+ik)) = PPQ% (4) 

■K~ ■= {'^"" iP2)\Ai,{A)\K+ {k)) = + Lf,R), (5) 

where P = pi + p 2 , Q = pi — P 2 , L = k — pi — p 2 - The form factors F, G, R and iJ are dimensionless scalar 
functions of the Mandelstam variables: 

s = {pi +P 2 f = {k- Lf, 

t = {k-pif = {P2 + Lf, (6) 

u = {k-p 2 )‘^ = {pi +Lf. 

We further define the invariant squared energy of the lepton pair Si = L^. For the hadronic matrix element, 
we regard Si as a fixed external quantity. 

2.2 Analytic Structure 

Let us first study the general properties of matrix elements of the hadronic axial vector current. It is instructive 
to draw a Mandelstam diagram for the process (see figures 1 and 2): since s + t + u = + 2M^ + si =: Sq 

is constant (for a fixed value of si), the Mandelstam variables can be represented in one plane, using the fact 
that the sum of distances of a point to the sides of an equilateral triangle is constant. 

The same amplitude describes four processes: 

• the decay K^{k) —>■ 7 r+(pi) 7 r“(p 2 )Ajj(L), 

• the s-channel scattering K'^{k)Af^{—L) —)• 7r+(pi)7r“(p 2 ), 

• the t-channel scattering Ar+(fc)7r“(—pi) —>• tt~ {p 2 )Aj^{L), 

• the u-channel scattering Ar+(fc)7r+(—P 2 ) —t 

The physical region of the decay starts at s = 4M^ and ends at s = {Mk — The s-channel scattering 

starts at s = {M^ + y^)^- If s^ = 0 is assumed, the two regions touch at s = (figure I). 

The sub-threshold region s < Sg := 4M^, t < to := {Mx + u < Uq := {Mx + forms a triangle 

in the Mandelstam plane where the amplitude is real. Branch cuts of the amplitude start at each threshold sg, 
to and Uq. There, physical intermediate states are possible (tttt intermediate states in the s-channel, Kir states 
in the t- and u-channel). 
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s-channel 



Figure 1: Mandelstam diagram for Ku for the case S£ = 0 
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s-channel 



Figure 2: Mandelstam diagram for Ku for the case S£ > 0 
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2.3 Isospin Decomposition 

Let us study the isospin properties of the Kn matrix element of the hadronic axial vector current in the different 
channels: we decompose the physical amplitude into amplitudes with definite isospin. 


2.3.1 s-Channel 

We consider the matrix element 


■K =(7r+(pi)7r {p2)\Af,{0)\K+{k)). 

As the weak current satisfies A/ = ^ the initial and final states can be decomposed as 


A^(0)|if+(fc)) = ^|l,0) + 4|0.0), 


r+Z 


V2' 


(7r+(pi)7r (p 2 )| = 


1 


1 


V3' 

1 


TT (Pi)7r+(P2) =^(2,0 -^(1,0 +^(0,0. 


V 6 






Hence, we can write the following decomposition of the matrix element into pure isospin amplitudes: 

A+- - -AA') + — 

2 ^ 76 ^ ’ 

2 ^ 76 ^ ■ 


Using AJ (fc, —L —)• pi,p 2 ) = A^~^{k, —L —> p 2 ,pi), we find the following relations: 


(7) 


( 8 ) 


(9) 


Al{k,-L ^ Pi,P2) = y^{A^ (fc,-L ^ pi,P2) + A+ {k,-L ^ p2,pi)), 
Al{k,-L ^ pi,P2) = (A+"(fc,-L ^-^ 1 ,^ 2 ) - A^~{k,-L ^ p2,pi)) . 


The pure isospin form factors are related to the physical ones by 


(s, t, u) = yi {F{s, t, u) + F(s, u,t)), 
(s, t, u) = yi (G(s, t, u) - G{s, u, t )), 

(s, t, u) = (i?(s, t, u) + i?(s, u,t)), 

= F{s,t,u) — F{s,u,t), 

(s, t, u) = G(s, t, u) + G(s, u, t), 

(s, t, u) = R{s, t, u) — R{s, u, t). 


We further note that 

-L pi,p2) = A^°\k, -L -)■ p2,pi), 
-L pi,P2) = -A^i^Hk, -L P2,Pi), 
and that the form factors of the pure isospin amplitudes satisfy 


(s, t, u) = (s, u, t), 

G^^\s,t,u) = —G^°^(s, M, t), 
R^^\s,t,u) = t), 


( 10 ) 


( 11 ) 


( 12 ) 
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FA'>(^s,t,u) = —FA'>(^s,u,t), 
G^^\s,t,u) = G^^\s,u,t), 
R^^\s,t,u) = —i?^^^(s, M, f). 


(13) 



2 . 3.2 t- and m-C hannel 


In the crossed t-channel, we are concerned with the matrix element 

A+~ = {n-{ p2)\Ai,{ 0)\K+{k)TT-{-pi)). 

In the u-channel, we analogously look at 

A+~ = {Tr+{pi)\Af,{0)\K+{k)TT+{-P2)). 


(14) 


(15) 


Note that due to crossing, these matrix elements are described by the same function - or its analytic continuation 
- as the corresponding s-channel matrix element. 

The t-channel initial and final states have the isospin decompositions 


|A'+(t),r-(-p.)) = \/||).-)) + /I 


■(P2) ^m(O) = \/ o ( 


2 /I 1 


3 \2 


3 1 

2 ’ 2 

1 /3 1 

3 \2’~2 


(16) 


(17) 


L,P2)- 


(18) 


whereas in the rt-channel, we are concerned with a pure isospin 3/2 scattering: 

\K+{k)TT+{-P 2 )) = 

(7r+(pi)|A^(0) = ^ 

We find the following isospin relation: 

■^^if^'^Hk,-p2 L,pi) = A^~ik,-L ^ pi,p2) 

= -)■ L,P2) + 

Note that the third component of the isospin does not alter the amplitude: just insert an isospin rotation matrix 
together with its inverse between in- and out-state to rotate the third component. 

The amplitude that describes pure isospin 1/2 scattering in the t-channel is then 

L,p 2) = ‘^A^^/‘^\k,-p2 L,pi) - ^A^i!^'^\k,-pi L,p2). (19) 

Defining analogous form factors for the isospin 1/2 amplitude, we find 

= ^F(s,f,u) - ^F{s,u,t), 

G(^/2)(s,t,M) = ^G{s,t,u) + ^G{s,u,t), (20) 

= ‘^R{s,t,u) - ^i?(s,'u,t). 

In the case = 0, it may be convenient to look at a certain linear combination of the form factors F and 
G, as we did in [12, 13, 14]: 


PL 

Fi := XF + {u- ^)^Gi 


( 21 ) 


where X := s,se), PL := \ — s — si) and A(a, b, c) := -I- 6^ -f — 2(a&-|- &c-|-ca) is the Kallen 

triangle function. 

Here, too, we can define the corresponding isospin 1/2 form factor: 

(s, t, u) := (s, t, u) + {u- ^) (s, t, u) 


PL 

= ( XF{s,t,u) + {u-t) — G{sA,u) 


1 


PL 


- o XF{s,u,t) + it-u) — Gis,u,t) 


2X 


( 22 ) 


3 1 

= 2^i{s,t,u) - -Fi{s,u,t)- 



2.4 Unitarity and Partial-Wave Expansion 

In this section, we will investigate the unitarity relations in the different channels and work out expansions of 
the form factors into partial waves with ‘nice’ properties with respect to unitarity and analyticity: the partial 
waves shall satisfy Watson’s final-state theorem. As we will need analytic continuations of the partial waves, 
we must also be careful not to introduce kinematic singularities. 

The derivation of the partial-wave expansion has been done for the s-channel in [16]. We now apply the 
same method to all channels. 


2.4.1 Helicity Amplitudes 

The quantities that have a simple expansion into partial waves are not the form factors but the helicity ampli¬ 
tudes of the 2 —)■ 2 scattering process [17]. However, helicity partial waves contain kinematic singularities. In 
order to determine them, we use the prescriptions of [18]. 

We obtain the helicity amplitudes by contracting the axial-vector-current matrix element with the polarisa¬ 
tion vectors of the off-shell W boson. In the W rest frame, the polarisation vectors are given by: 


£^ = ( 1 , 0 , 0 , 0 ), 

4 = ^(0,0,±l,t), (23) 

4 = ( 0 , 1 , 0 , 0 ). 


Si = 


/ 0 0 0 
0 0 0 
0 0 0 
\ 0 0 I 


0 \ 

0 

—i 

0 J 


= sf 


^ and Si, defined by 


/ 0 0 

0 0 ' 


0 0 

0 i 


0 0 

0 0 


\0 -i < 

0 0 


1 

/ 0 ( 

+ 4' + si = 

0 : 
0 ( 


1 

^0 ( 

are 

: listed below: 


A A 

^0 

s 

0 1 

1 

Si 

0 ±1 

0 


S3 = 


0 0 0 ^ 
2 0 0 
2 0 
0 0 2 / 


/ 0 0 
0 0 
0 i 
\ 0 0 


0 

—i 

0 

0 


0 
0 
oy 


(24) 


If we boost the polarisation vectors into the frame where the W momentum is given by L = (L°,L^,0,0), 
= Si, we obtain: 


e(‘ = ^(L°,Li,0,0), 

l/Si 

£l = T(o,o,±l,!), (25) 

The contractions of these basis vectors with give the different helicity amplitudes: 

A := Aef • (26) 

We extract the kinematic singularities by applying the recipe of [18], chapter 7.3.5, to these helicity ampli¬ 
tudes. 
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2.4.2 Partial-Wave Unitarity in the s-Channel 
2.4.2.1 Helicity Partial Waves 

The unitarity relation for the axial vector current matrix element reads 

lra(iA^P{k,-L ^ pi,P2)\ [ dqidq2{2TT)‘^6^'^\pi + P2 - qi - q2) 

^ ^ 4 J (27) 

T^^^*(qi,q 2 -)■ Pi,P 2 )iA[^\k,-L ^ < 71 , 92 ), 

where dq := (■ 2 ,f) 32 g 0 is the Lorentz-invariant measure and where a symmetry factor 1/2 for the pions is included. 

denotes the elastic isospin I Trir-scattering amplitude. Note that this relation is valid in the physical region 
and that kinematic singularities have to be removed before an analytic continuation. 

We perform the integrals: 

lm(^iA[^\k,-L ^ pi,P 2 )^ = J d^” {s,cose')iA^P {s,cosO'' (28) 

where a^ris) = ■\/l — s and of course cos 9' has to be understood as a function of cos 6” and (j)" through 
the relation 


cos 6' = sin 9 sin 9" cos (j)” + cos 9 cos 9". 


(29) 


If we expand T and Ai into appropriate partial waves, we can perform the remaining angular integrals and find 
the unitarity relations for the it '^4 partial waves. 

We expand the Trir-scattering matrix element in the usual way: 


(s, cos 9') = ^ Pi (cos 9') tj (s) 




with 


tl{s)= tlis) 


(30) 


(31) 


The Ki 4 ^ helicity amplitudes are expanded into appropriate Wigner d-functions, which satisfy dgQ (9) = 
Pi(cos9) and <iig(d) = —[l{l + 1)]“^/^ sindP/(cosd). We have to take care of the kinematic singularities of the 
helicity amplitudes [17, 18]: 


iAY\s, cos 9) =J2 Pi{cos9) j ^2 


(s)o'7r(s)\ (/) 

“t.i 


(s), 


1=0 


iA[^\s, cos 9) = iAY’ ''^^^' = 


; aP) (^) _ ^Kl 


K 

1/2/ 00 


z.42^^ (s, cos 9,4)) = iA 2 ^^ sin 0 = sin d ^ Pi (cos 9) 

1=1 


- aP) . 


J2Piicos9) 
1=0 

,1/2 


^Ki ^ (I) 


Ml 


(^)<^'n' (^) 


/-I 




COS 4 0,^1 (s), 


Ml 


(32) 


where A^'^ := .4^^ — .4^^ The square roots of the Kallen function cancel exactly the square root branch cuts 
in the Legendre polynomials between [Mk — and {Mk + y^)^- factor Ml in the denominators 

appears only for dimensional reasons. All the defined partial waves afi are free of kinematic singularities and 
can be used for an analytic continuation from the decay region through the unphysical to the scattering region. 

If we insert the partial-wave expansions into the unitarity relation, the remaining angular integrals can be 
performed and the unitarity relation for the Kn partial waves emerges (i = t, 0, 2): 

Im (a|y(s)) = (33) 

In particular, we find that the phases of the K (^4 s-channel partial waves are given by the elastic TTTr-scattering 
phases (this is Watson’s theorem) for all s between 4Mj and some inelastic threshold: 






Ml (s 


(34) 
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2.4.2.2 Partial-Wave Expansion of the Form Factors in the s-Channel 

In order to find the partial-wave expansions of the form factors, we write explicitly the helicity amplitudes 
(generalised to a generic </>): 

1 /I, 


MKy/^i - S - + 2'^7r{s)X][^{s) cos6» -h se 


i 

Cq 




(35) 


A AR _ iAA)p^^ _ i AiApt^ 


V^CT,r (s) sin 6* cos (j) g^-^M. 


-72 

Since the contribution of the form factor R to the decay rate is suppressed by it is invisible in the 
electron mode and we do not have any data on it. We therefore look only for linear combinations of the form 
factors F and G that possess a simple partial-wave expansion. We find: 

pH) p <J.{s)PL{s) ^ ^ 

X(s) XKt{s) 


“S<»- 


1=0 
Mk 


(36) 




Z-1 


4?(s)- 


We write the partial-wave expansions of F and G in the form: 


r<'> = f; nicos 9) ( ' f,n (,) _ 9 g<'>, 

(=0 \ ^ / 


G''^=E7'(cos7 




(37) 




1^1 


Ml 


where the partial waves and g\^'^ satisfy Watson’s theorem in the region s > AM"^: 


f\s)= f\s) 


piSi [a) 


9?\s)= g^l\s) 


jA («) 


(38) 


2.4.3 Partial-Wave Unitarity in the t-Channel 


2.4.3.1 Helicity Partial Waves 

The discussion in the crossed channels is a bit simpler because we are interested in partial-wave expansions 
only in the region t > {Mk + Mj^Y or u > {Mk + MAl, i.e. above all initial and final state thresholds and 
pseudo-thresholds. Therefore, we do not have to worry about kinematic singularities, since we will not perform 
analytic continuations into the critical regions. 

In the crossed channels, we consider Ktt intermediate states in the unitarity relation: 


(39) 


Im ^ L,P 2 )) J dqKdq^{2Tr)'^d^^\k-pi-qK-q,T) 

{qK,qiA k,-pi)iA[^^'^\qK,q7T L,P2), 

where jg the isospin 1/2 elastic TfTr-scattering amplitude. By performing the integrals we obtain: 

lm(iA[^^^\k,-pi T,P 2 )) = J cos e{)iA[^^^\t, cos (40) 
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The Kt: scattering matrix element is expanded in the usual way: 

OO 

T^^^‘^^t,cos9t) = y^P;(cos 

1=0 


with 


= iy\t) .«;'■<■> 


I 

We expand the Ku helicity amplitudes as follows: 

= f;fl(cos9,) 

1=0 \ ^ J 

OO 


(41) 

(42) 


(43) 


= ±mcose,) f)(t), 

1=0 \ ^ / 

cos 6t,(j)t) = iA[^^‘^\t, cos 61 , (j)t) — iA^y^\t, cos 9t, (j)t) 

By inserting these expansions into the unitarity relation (40), we find that all the partial waves satisfy Watson’s 
theorem {i = t, 0, 2): 


<•!/’(*) = 


-Tm 




(44) 


2.4.3.2 Partial-Wave Expansion of the Form Factors in the t-Channel 

By contracting the axial vector current matrix element in the t-channel with the polarisation vectors, we find 
the helicity amplitudes (for a generic As we are not interested in R, we do not need the component: 

(A^^(t)(M2 -si-t) cos9t + (<)(M|- - +1)^ 


K^Sj; y4t \ “ / 

+ ^ -se- t)cos9t + A]/^(t)(M|. - - 3t)) ^ 

/I/O \ 


(45) 


= tA<i^/^^el - tAl!/^h^_ = f sin0, cosc/>t + G''/^)) 

y2Mx \ \/t ^ ^ 

This results in the following partial-wave expansions of the form factors: 

riv^> = f;a(c««,) 


1=0 

1 

n 


aK(() 


{Ml -Si- cos 9t j ^ F/ (cos^*) 

(G / 1=1 




M% 


G^y-)=-Y^R{cos9i) /44t) 


(46) 


+ l^{Ml-Ml+t+{Ml 


g ) J 1=1 \ ^ / 


where also the new partial waves and satisfy Watson’s theorem in the region t > {Mk + AG)^' 

n/2'l/.N I ('I/2 'I..nI 




(47) 
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2.4.4 Partial-Wave Unitarity in the i6-Channel 


2.4.4.1 Helicity Partial Waves 

The u-channel (i.e. the isospin 3/2 case) can be treated in complete analogy to the t-channel. We start with 
the unitarity relation: 


Im ^ J dqKdq^{27r)‘^S^'^\k - P 2 - qK - q^:) 

■ {qK,q-K k, -P 2 ) iA[^^‘^\qK, q^ ^ T,pi) 


1 1 '^if7r(^) 


8 (2^)2 

The it'TT-scattering matrix element is expanded as 


J dn'^T^^^^'>*{u,cos0'JiA[^^‘^\u,cose”,(j>'^ 


(48) 


T^^^‘^\u,COS0u) = y^P;(cOS 0u)tf'^{u), 
1=0 

P/2^„A _ L3/2. 


tr\u) = tr\u) 


3 ./ 2 ( 

pi^V^A) 


and the itr ^4 helicity amplitudes according to 


cos 0u) 

00 / 

= ^Pl{cOS0u) 1 
1=0 ' 


Mf, 

iAl^^‘^\u, cos 0u) 

00 / 

= '^Pl{cOS0u) ( 
/=0 ' 


< 

! 2 ^‘^\u, cos 0u,(t>u) 

= /.4^|i^^^^(u,cos0, 



<r{u\ 


sin 9. cos.^. f; p;(co. 0 .) ( 


Z-1 








(49) 


(50) 


Performing the angular integrals in the unitarity relation, we find that the partial waves satisfy Watson’s theorem 

(* = f,0,2): 


Im 




1 1 A 


1/2/ 

kA 


2/ + 1 IOtt u 


,3/2*, X (3/2), S 
h AAiJ '{u), 


J3/2) 

^i,l 


(u) = 


,(3/2) 


(u) 




2 . 4 . 4.2 Partial-Wave Expansion of the Form Factors in the u-Channel 

The contraction of the axial vector current matrix element with the polarisation vectors yields: 


*-4o^/2) = = Mk^i, (^ 4 /) -Si-u) cos6i„ + A]/^(u)(M|, - + m)) F 

- ^ -Si-u) COS0U + A]/^(m)(M|. - - 3m)) , 

£^_ = ^ sin6l„ cos(/„ (F — G) V 
V2Mk Y yii J 


1 


(51) 


(52) 
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Hence, the partial-wave expansion of the form factors is given by 


F = ^Pi{cos0u) 




1=0 


Mi 


- ^(mk - M^-3u+ {M^- Si- cos6»„ 




G = y^Fjcos 6u) 


1 .(3/2) 


(53) 


(=0 

1 I 

2u ' 


Mi 




Ml-Ml +U+ (Ml -Si-u) cos f; P/(cos ^(3/2) 

^ ^TT W 


Z=1 


Mi 


where the partial waves and satisfy Watson’s theorem in the region u > {M^ + 


, gr^>iu) = \gr^>iu) 

2.4.5 Projection and Analytic Structure of the Partial Waves 

The several partial waves and g^^ can be calculated by angular projections: 

'"'HAf .. 

where z) := xG')(^s^ t{s, z), u(s, z)), X € {F, G}, I G {0,1} and 

t{s, z) = ^ (Eo - s - 2Acr^z), 
u{s, z) = ^ (Eo - s -I- 2Acr^z). 

Since t{s, —z) = u{s, z), the definition of the pure isospin form factors (11) implies 

fi°^ (s) = 5^°^ (s) = 0 y I odd, 

=0 y I even. 

Hence, we can as well directly use the partial waves of the physical form factors: 

f Ml y 21 + 1 f\ , , <jAs)PL{s) 


p(3/2). 




(3/2), 


J3/2), 


i6^ (li) 


Ml 

gi{s) = ' 


1-1 


,XKi{s)a^{s)J 
which still fulfil Watson’s theorem 

Ms) = \Ms)\F^i^^^ , gi{s) = \gi{s)\F^!(^\ 


f’i-i(^) Pi+iiz) 


(54) 


(55) 


(56) 


(57) 


(58) 


(59) 


where I = {I mod 2). 
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In the crossed channels, the partial wave projections are given by 




K 


+ ^ ~ — i + — Si — j 


21 + 1 




'-1 


V 




(1/2)/^n _ [ _ Mk _ 


1-1 




[' ^^^Pi-ii+)-Pi+ii+) F^^/^\t,zt) + G^^/^\t,zt) 
-1 


^(3/2)/ X 

fi m = 


Mf, 


21 + 1 




dZuPl{Zu) 


2 2 
F(u,Zu) + G{u,Zu) 


1-1 


+ L(mI,-M^-u+ {Ml - Si -u) ~ ^ 


2u 

\ 

(3/2)/ N 
9i {u) = 




Mj, 


1-1 




Pl-l {Zu) - Pl+l{Zu) F{u, Zu) - G{u, Zu) 


1-1 


2 2 

where XP'>{t,Zt) := XP'> {s{t, Zt),t,u{t, Zt)), XP\u,Zu) := XP'> {s{u, Zu),t{u, Zu),u), X € {F,G} and 


(60) 



The construction of the partial waves has been done in a way that excludes kinematic singularities for 
s > 4:Ml and t,u > {Mk + There may still be kinematic singularities present below these regions, but 

they do not bother us. But also the analytic structure of the partial waves with respect to dynamic singularities 
is not trivial. 

For the s-channel partial waves, there is of course the right-hand cut at s > AMl. Further cuts can appear 
through the angular integration, i.e. at points where the integration contour in the t- or u-plane touches the 
crossed channel cuts. If s lies in the physical decay region, the integration path is just a horizontal line from one 
end of the decay region to the other (see the Mandelstam diagram in figure 2). When we continue analytically 
into the region {Mk — y/se)'^ < s < {Mk + v^)^> ^^e integration path moves into the complex t- and u-plane 
and crosses the real Mandelstam plane &t t = u: the square root of the Kallen function X = \Xyi{s) is purely 
imaginary in this region. One has to know which branch of the square root should be taken. The correct sign 
is found by taking s real and shifting Mk —> Mk + ie (see [19]). With this prescription, the Kallen function 
turns counterclockwise around Xki = 0 when s runs from s < {Mk — y^)^ to s > {Mk + y/si)'^- The square 
root of the Kallen function therefore takes the following values: 

f +|A^/"(s)| s<{MK-V^if, 

A^^'(s) = < +^|A^/"(s)| {Mk - ^i)^ < s < {Mk + (62) 

[ -|A^/,"(s)| {MK + V^i)^<s. 


In the region s > {Mk + y/si)'^, the integration path again lies in the real Mandelstam plane from one to the 
other end of the scattering region. 

As we are away from the t- and u-channel unitarity cuts at t,u > {Mk + Mn-)^, this extension of the 
integration path into the complex t- and u-plane is the only subtlety that has to be taken into account. 

In the region s < ‘IMl, there is a left-hand cut at s € (—oo,0): the integration path extends again into the 
complex t- and u-plane in the region 0 < s < AMl (due to the second square root). It diverges at s = 0 and 


15 



returns to the real axis at s < 0, but this time it touches the t- and w-channel unitarity cuts at t,u > 
which produces the left-hand cut of the s-channel partial waves. 

This left-hand cut can be most easily found by looking at the end-points of the integration paths: solving 
the equation 


= Mzp = ^ (So - s=F 2X(s)(j^(s)) (63) 

for t± > {Mx + gives the left-hand cut s € (— 00 , 0). 

Let us consider the crossed t-channel (the situation in the it-channel is analogous). We have defined the 
partial-wave expansion in the scattering region t > {Mk + Tf^)^. Therefore, we also define the square root 
branches of the Kallen functions and in this region. The sign of the square root branch can be absorbed 
into the definition of the partial waves. 

The right-hand t-channel unitarity cut at t > {Mx + also shows up in the partial waves. A second 

possibility for singularities in the t-channel partial waves arises when the integration path touches the s- or 
M-channel unitarity cuts. For t > {Mx + Mn-)^, the integration path lies on the negative real axis of the 
s- and M-planes (this can be seen in the Mandelstam diagram in figure 2). In the region {Mx — Mn-)^ < 
t < {Mx + the integration path extends into the complex s- and w-plane. For the value of t fulfilling 

h (Eq — t — jAxTr^iTr) = 4M^, the integration path in the s-plane touches the s-channel branch cut. From this 
point on towards smaller values of t, the integration path has to be deformed in the s-plane. Since the u-channel 
cut appears only at u > {Mx + Mn-)^, such a deformation is not needed in the u-plane. At t = {Mx — Mt^Y, 
the integration path in the s-plane has then the shape of a horseshoe wrapped around the s-channel cut. For 
even smaller values of t, the path unwraps itself in a continuous way, such that for t < \{Mx — ‘2,M‘^ + si), the 
integration path lies completely on the upper side of the s-channel cut. 

The cut structure in the t-channel partial wave is rather complicated, at least for si > 0: The left-hand cuts 
can be found by solving the equations 

■s± = X f^o - i + - (± A]/^(t)Aj)^^(t) - , 

t ^ (± - Aic.A,^) j , 

for s± > and u± > {Mx + MtH. While the second equation results in a cut along the real axis, the first 
equation produces an egg-shaped cut structure in the complex t-plane with Re(t) < {Mx — M^^)'^, shown in 
figure 3. The exact shape depends on the value of s^. 



Figure 3: The left-hand cut of the t-channel partial waves {se = 0.3M^). 
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2.4.6 Simplifications for 0 


In the experiment, a dependence on has been observed only in the first partial wave of the form factor F 
[6, 9]. If we neglect this dependence on si and assume that = 0, the treatment can be significantly simplified. 


• The square root of the Kallen function simplifies to 

= Afi - s, 

S £—>-0 

the square root branch cut disappears. Hence, the integration path for the angular integrals in the 
s-channel always lies on the real axis. 

• The left-hand cut structure of t- and u-channel partial waves simplifies to a straight line along the real 
axis. The egg-shaped branch cuts disappear in the limit Si —>■ 0. 

• From (37), we see that the quantity 


lim Fy'^ = lim I - 

st->-0 Sf-J-O ' 




Mk^pd) 


u — t 


g(f 


(65) 


has a simple s-channel partial-wave expansion into Legendre polynomials. If we consider (46) in the limit 
Si —)■ 0, we find that exactly the same linear combination of the form factors and has a 

simple t-channel partial-wave expansion into Legendre polynomials. The same follows from (53) for the 
u-channel. In this limit, the form factor Fi can therefore be treated independently from the other form 
factors. This is the procedure that has been followed in [12, 13, 14]. 


There are several reasons why we abstain here from taking the limit Si —)■ 0, which would result in a 
substantial simplification of the whole treatment. The experiments provide data on both form factors F and G. 
In order to include all the available information, we deal with both form factors at the same time. There is also 
some data available on the dependence on si, which we include in this treatment. And finally, the matching to 
yPT becomes much cleaner if it is performed with F and G directly, since these are the form factors with the 
simplest chiral representation. 


2.5 Reconstruction Theorem 

Since the form factors F and G describe a hadronic four-‘particle’ process, they depend on the three Mandelstam 
variables s, t and u and therefore possess a rather complicated analytic structure. However, it is possible to 
decompose the form factors into a sum of functions that depend only on a single Mandelstam variable - a 
procedure known under the name of ‘reconstruction theorem’ [20, 21[. Such a decomposition provides a major 
simplification of the problem and leads to a powerful dispersive description. 


2.5.1 Decomposition of the Form Factors 

The explicit derivation of the decomposition of the form factors F and G into functions of a single Mandelstam 
variable can be found in [15[. It is based on fixed-s, fixed-t and fixed-M dispersion relations. We have to assume 
a certain asymptotic behaviour of the form factors, e.g. for fixed u, we assume 


lim 

|s|—>-oo 


s" 


lim 

|—>-oo 


xrit) 


= 0 , 


( 66 ) 


where the Froissart bound [22[ suggests n = 2. However, we are also interested in the case n = 3 in order to 
meet the asymptotic behaviour of the NNLO xPT form factors. We therefore write down either a twice- or 
thrice-subtracted dispersion relation for the form factors. Then, we use the partial-wave expansions derived in 
the previous section. We neglect the imaginary parts of D- and higher waves, an approximation that is violated 
only at 0{p^) in the chiral power counting. It implements the case sg yf 0. 
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The result of the decomposition is the following: 


F{s,t,u) 


G{s, t, u) 


Mo{s) + FIi{s) 


3 ^ ^ 3 ^ ^ 3 2M| 

1 n /i\ ^ ~f ^K'k^I'K 

+ -774- 


iVi(t) 


Mi 




I p / \ I 0 “f ^-fCTT^^TT p / \ ^KtT 3m p / n 

+ ^o(m) + -^^^4-- 9A72. 


2M|, 


+ o(/), 

Mi(s) 


7^oW - 7-- Ni[t) + 7 ^^^2 


3 3 

1 


3 2M|. 


+ -RoIm) + - Jj-4^ - Rl{u) - 7772 — i?i(u) 

+ 0(/)- 


2M|. 


In the case n = 2, the various functions of one variable are given by 


Mo(s) 

Mi(s) 

Mi(s) 

No{t) 

N,{t) 

iVi(t) 

Ro{t) 

Ri{t) 

Ri{t) 





Im/o(s') 

(s' — s — ie)s'^ 


ds', 


TOi H- 


rr^-^.4/.(.')-?^^».(¥)^.'. 

Ac (s'-s-ie)s' V Aa:£(s') 7 


+ m\ 




sq 

s 

Ml 

!! 


{s' — s — ie 

TT 


roo 


Iso {s'-s-ie)s 

lm/f'¥0 ,,, 


/2 


ds', 


1 

TT 


*0 


f — t — ie 


to it' - t- ie)t'‘^ 


t /•- M| 

TT t' [f -t- ie)t' 



Im/f/^^(t') 

{f — t — ie)t'^ 


dt', 


1 r ‘ Im 

TT 7to t' -t- le 


(/f'">(i')+ 


(^r-TT + t')Fll 

2t'XeAt') 



t_ Ml lTngf/^\t') 

TT t' [f -t- ie)t' 


(67) 


( 68 ) 
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while for n = 3, the functions of one variable are 


Mo{s) 

Mi(s) 

Mils) 

Noit) 

Ni{t) 

Niit) 

Ro{t) 

Riit) 

Riit) 


0 , 1 S , 2 /■“ Wo(s') J / 


n ^ ^ 


So 

^2 


{s' — s — ie)i 


(s' — s — ie)s 


^Kt{s') 




+ TOi 

t3 


°M2 


, s3 /■“ Im 5 i(s') ^ , 
(s'-s-ze)s'3^"’ 

■dt', 


,so (s'-s-ze)s 


^ iio {f -t- 

^ + - r ^ Im f + ( 1 / 2 ). /A ^ , 






'to 

(3/2), 


t' {f — t — ie)t' 




dt', 


11 Jto {f -t - ie)t'' 
t f°° 1 

71 Jto {f -t- ie)t‘ 

TT Jto t' {t'-t-te)t'^ ■ 




(69) 


Actually, since the P-wave of isospin / = 3/2 A'tt scattering is real at 0{p^), so are the partial waves 
and Hence, the functions Ri{t) and Ri(t) could be dropped altogether in the decomposition. The phase 

3 /2 

61 is also known to be tiny in phenomenology. 


2.5.2 Ambiguity of the Decomposition 

We have decomposed the form factors F and G into functions of one variable Mo(s), .... However, while the 
form factors are observable quantities, these functions are not. It is possible to redefine the functions Mo(s), 
... without changing the form factors and hence without changing the physics. 

Therefore, let us study this ambiguity of the decomposition of the form factors. We require the form factors 
to be invariant under a change of the functions of one variable: 

Mo{s) I—)■ Mq{s) + 5Mo(s), 

Mi(s) Mi(s) + 6Mi{s), (70) 
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which we call ‘gauge transformation’. The shifts have to satisfy 


0 = 6Mo{s) 


u — t 


SMi{s) 


+ IsNoit) + 


K 


3 2M 


K 


6Ro{u) 4 


1 f(s m) + 1 Ai^-.^. 3t ^ ^ 

‘ 3-UJ- 

u{s “1“ c r> / \ ^K-k 3^/- c / \ 

-^- SRi{u) -— SRi{u), 


M 


K 


K 


0 = SM^{s) 

- - 3 


2 t[s u) + Ax-JT^i-JT f; 2 Axtt t XAT u\ 

- Mi -+ 3^sr ’ 

f4.\ ^ tis ~ u) Ax-jrAij^ lAxTT+t ^ ^ ^ 

- - 3- Ml - 


(71) 


(72) 


\ XT? / \ \ t) XT? ( \ U ~ . , 

4 i5i?o(w) +-- 6Ri(u) - 2M^ 5Ri{u). 

The solution to these equations is found in the following way: we substitute one of the three kinematic variables 
by means of s + t + u = Eq- Then, we take the derivative with respect to one of the two remaining variables 
and substitute back Eq = s + f 4 u. After four or five such differentiations, one gets differential equations for 
single functions 5Mq, ... with the following solution: 


SMois) = 

II 

Tcfosd 

Mn 2 

- C 2 S 

? 


SMi{s) = 

_ Mt 
- ^0 

4cf’s4 

Ml 2 

- C 2 S 

? 


6Mi{s) = 

_ 

1 Ml , 

-f ^5 4 


4 

Ml 3 
C 3 's 

SNoit) -- 

= 

-1 , ^No 
-h Cq 

4 cf»t 

4 


SNi{t) -- 

= c^\t 

-h Cq 

4 cfH 



6Ni{t) = 

= S,t 

-h Cq 

4 cfH 

4 

c§Hl 

SRo{t) = 

= 

-1 1 -Ro 
-h Cq 

+ cf“t 

+ 


SRi{t) = 

= C^\t 

-h Cq 

+ 0^4 



SRi{t) = 

= C^{t 

-h Cq 

+ c^H 

+ 



(73) 


Inserting these solutions into the various differential equations results in algebraic equations for the diverse 
coefficients. In the end, there remain 13 independent parameters. In complete generality, we therefore have a 
gauge freedom of 13 parameters in the decomposition (67). The gauge can be fixed by imposing constraints on 
the Taylor expansion or the asymptotic behaviour of the functions Mq{s), .... 

First, let us restrict the gauge freedom by imposing the same vanishing Taylor coefficients as in (68), i.e. we 
exclude all the pole terms, the constants in Nq, Ni, Rq, Ri and even a linear term in Rq. Then, we further 
demand that asymptotically the functions behave at most as in (69), i.e. like Mi{s) = 0(s), Mi(s) = O(s^), 
Ni(t) = 0(1), Ni{t) = 0{t), Ri(t) = 0(1) and Ri{t) = 0{t). After imposing these constraints, we are left with 
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a restricted gauge freedom of 3 parameters, which we call and 

SMois) = (2A^^ - + 2C7^«) ~ "" ’ 

m,(s) = - (a’- + B”- + 2C«-) ^ + B*- Ig; + (b*. + 2C«-) 

Sm{s) = Eg _ (2A^i + _ 26-^0) Ak.A,^ + B^^^^Ak.) ^ 

2 

- + B^^ - 2C7««) 2Eo) + (2^^^ + B^^ - 2C^») ^ 

5Ao(t) = - (2^«i - + 2C'««) " 10C'««) 

,5Ai(t) = - i (2A^i + 35^1 - 6(7^“) , 

5iVi(t) = - (6^«i + + 6C'^“) 

SRi(t) = A^\ 


(74) 


In order to fix the gauge completely, we have to impose further conditions. We will use two different gauges. 
The first one corresponds to the case of an asymptotic behaviour that needs n = 2 subtractions. It is most 
suitable for our numerical dispersive representation of the form factors and for the NLO chiral result. In this 
case, the asymptotic behaviour excludes quadratic terms in 6Mq(s) and 6Mi(s) or a linear term in SMi{s). 
Hence, in the representation (68), the gauge is completely fixed. 

The chiral representation, being an expansion in the masses and momenta, does not necessarily reproduce 
the correct asymptotic behaviour. The 0(p®) chiral expressions show an asymptotic behaviour that needs 
n = 3 subtractions. In this case, one has to fix the gauge rather with the Taylor coefficients, e.g. by excluding 
a quadratic term in i? 0 i a constant term in i?i and a linear term in Ri. Therefore, also in the representation 
(69), the gauge is completely fixed. 

Note that the second representation (69) makes less restrictive assumptions about the asymptotic behaviour. 
Therefore, the first representation (68) is a special case of the second (69). One can easily switch from the first to 
the second representation with the help of the gauge transformation (74). In this case, the additional subtraction 
constants will be given by sum rules. 


2.5.3 Simplifications for 0 

As a test of the decomposition, let us study the linear combination 

Fi{s,t,u) = - s)F{s,t,u) + ^{u-t)G{s,t,u) 

in the limit —>■ 0. We neglect the contribution of the isospin 3/2 P-wave: 


(75) 




u — t 


lim Fi(s, t, u) = lim —A- F(s,t,u)-\ - G(s,t,u 

Mi - s . , 

= ^^2 - Mois) 


+ 


u — t 




+ g [(^ “ ^n)^oit)] + 2 [(^ “ Fl!^)Fo(i)] + (u- MI)Rq{u) 






•^K 


(76) 
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By identifying 


rFw„^ “ s 


M^^is) = 




Mo{s), 




N^^{t) = {t-M^)No{t), 
R^-{t) = {t-Ml)R^{t), 


(77) 


2M2. 


Ni{t), 


we recover the decomposition of the form factor used in [12, 13, 14]. We further note that the imaginary 
parts of the functions of one variable in this decomposition are given by 


ImM,f^ (s) = 


Mk-s 


Im/o(s), 


ImMfi(s) = 

Im<^(t) = (t-M2)Im/(i/2)(t), 

= {t- 

it) = 


(78) 


Tm N[^ 




and repeat the observation of section 2.4.6 that in the limit si —)■ 0, these partial waves are given by projections 
of Fi in all three channels. Hence, the form factor Fi decouples in this limit and can be treated independently 
in the above decomposition. 


2.6 Integral Equations 

2.6.1 Omnes Representation 

The decomposition of the form factors (67) signifies a major simplification, since we only have to deal with 
functions of a single Mandelstam variable. These functions (68, 69) are constructed in such a way that they 
only contain the right-hand cut of the corresponding partial wave. Their imaginary part on the upper rim of 
their cut is given by 


lmiMo(sj = lm/o(s), 


ImMi(s) = Im (^/i(s) - g^[s)^ , 

ImMi(s) = Im(/i(s), 

ImiVo(t) = Im/^^/^^(t), 

ImJV.(() = Im 

ImJV,(() = Im(^^9l''"’(')). 

ImRo(0 = Im/o^^^V). 

ImS.d) = Im 

Im«.(() = Im(^Misr>(i)). 


(79) 
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Therefore, we can write 


Mo(s) + Mo{s) — fo{s), 

M,{s) + Mi(s) = h{s) - 

>^Ke{s) 

Mi{s) + Mi(s) = 5 i(s), 

Ni{t) + Ni{t) = 




(80) 


Ml 


R,(t)+ s.(t)= 


2t\Ut) 


Ri{t) + ^i(t) = 


where the ‘hat functions’ Mq{s), ... are real on the cut: indeed, they do not possess a right-hand cut, but 
contain the (possibly complicated) left-hand cut structure of the partial waves (see section 2.4.5). Writing 
Im/o(s) = /o(s)e“*‘^o^®^ sin^J^s), ... leads directly to the following equations: 

ImMo(s) = (Mo(s) -I- Mo(s))e“*‘^®^®^ sin(5o(s), 

ImMi(s) = {Mi{s) + Mi(s))e“*‘^i^®^ sin(5i(s), 

ImMi(s) = {Mi{s) -|- Mi(s))e“*‘^i^®^ sin(5j(s). 


(81) 



it) 

• Jvl/2 

Sin 

(t), 


(t) 

• Jvl/2 

(t), 


it) 

• Jvl/2 

smdi 

(t), 


it) 

• j:3/2 

sm dg 

it), 


it) 

• s:3/2 

smdi 

it), 


it) 

• j:3/2 

it). 


where, below some inelastic threshold, the phases 6j agree with the elastic tttt- or iilTr-scattering phase shifts. 
Therefore, the functions Mq, ... are given by the solution to the inhomogeneous Omnes problem. The minimal 
number of subtractions appearing in the Omnes representation is determined by the asymptotic behaviour of 
the functions Mq, ... and the phases (5/. 

Let us extend these equations even to the region above inelastic thresholds by replacing 5 i—>■ w. 


ImMo(s) = (Mo(s) -I- Mo(s))e sina;o(s). 


where ujj (s) = 6j (s) -I- ry/(s) and gj (s) = 0 below the inelastic threshold s = h?. 
We define the usual once-subtracted Omnes function 


r2(s) := exp 




{s' — s — ie)s' 



(82) 


(83) 


If the asymptotic behaviour of the phase is lim 5{s) = rmr, the Omnes function behaves asymptotically as 

S—¥(X) 

0{s~'"'). Provided that the function M{s) behaves asymptotically as (P(s^), we can write a dispersion relation 
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for M/Q, that leads to a modified Omnes solution 


M(s) = n(s) 


Pn-l{s) H- 

TT 


M{s') sin(5(s') 
|n(s')|(s' — s — ie)s' 


rds' 



Mis') sin(5(s') , , 

- fi 

\n{s')\{s'-s-ie)s"' 

(M(s') + ReM(s')) sin 77 ( 5 ') 

- (I S 

|n(s')| cos((5(s') + ? 7 (s'))(s' — s — ie)s"'' 


(84) 


where the order of the subtraction polynomial is n — 1 = k + m. 

Actually, we do not know the phase S at high energies. Inelasticities due to multi-Goldstone boson interme¬ 
diate states, i.e. more than two Goldstone bosons, appear only at 0{p^) [20], hence the most important inelastic 
contribution would certainly be a if AT intermediate state in the s-channel. This could be included by using 
experimental input on 77 up to s ss (1.4 GeV)^. 

We could make a Taylor expansion of the inelasticity integral and neglect terms that only contribute at 0{p^) 
to the form factors by applying the power counting ^ P^- This would introduce quite a lot of unknown 
Taylor coefficients. Here, we follow another strategy: we set 77 = 0 and assign a large error to the phases 5 at 
high energies. We assume further that 6 reaches a multiple of tt above a certain s = A^. The two high-energy 
integrals in (84) drop in this case. 

Assuming that the phases behave asymptotically like dg — )■ tt , (5( — > tt and all other Sj —> 0, we find the 
following solution for the case of n = 2 subtractions: 


Mo{s) 

Mi{s) 

Mi{s) 

No{t) 

Niit) 

Niit) 

Roit) 

Ri{t) 

Ri{t) 


n}(s) 

n}(s) 

^Tit) 

^Tit) 


TT 


Mi 


Mi 


Mo(s')sin(5g(s') 
|f 2 Q(s')|(s' — s — ie)s'^ 


ds' 


,Mi 


. 5 M 1 ^ ^ ^ 
Ml ^ TT , 


Ml (s') sin i5((s') 
|ni[(s')|(s' — s — ie)s'‘ 


:ds' 


Ml (s') sin i5i(s') 
|f2i[(s')|(s' — s — ze)s'^ 

p t'\, 

’’’-'io ^{t')\{t'— t — ie)t'^ / 

1 Ni{t')smSy‘^{t') 

^ Jto — t — ie) I 


Ml 7 ,Mi ^ ^ 

“ + Mi + Ml ^ TT , 



^ Jto — t — ie)t' I 

^ Rojt') sin it') 

^ ./to |ng^^(t')|(t' — f — 7e)f'^ i 

1 Ri{t')sinSl^‘^{t') 

^ Jto |n5^^(t')|(t'— t — ze) J 


^ Jto — t — ie)t' i 


(85) 


where we have fixed some of the subtraction constants in Nq, Ni, Rq and Ri to zero by imposing the same 
Taylor expansion as in the defining equation ( 68 ). 

Note that driving the Ktt phases to 0 is somehow artificial. They are rather supposed to reach tt at high 
energies. However, this would introduce three more subtraction constants in our framework. Since the high- 
energy behaviour of the phases does not have an important influence on our results, we abstain from introducing 
more subtractions and take these effects into account in the systematic uncertainty. 

In the case of n = 3 subtractions, six additional subtraction constants appear in the Omnes representation. 
The conversion of a solution for n = 2 into a solution for zz = 3 requires a gauge transformation in the Omnes 
representation, as explained in appendix G.l. 
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2.6.2 Hat Functions 


The hat functions appearing in the Omnes solution to the functions of one variable (85) can be computed 
through partial-wave projections of the form factors: (80) should be understood as the defining equation of the 
hat functions. One has to compute the partial-wave projections of the decomposed form factors F and G (67) 
and subtract the function of one variable (Mg, ...). Finally, one obtains an expression for the hat functions 
in terms of angular averages of the single-variable functions. The explicit expressions for the hat functions are 
given in appendix C.2. 


3 Numerical Solution of the Dispersion Relation 

3.1 Iterative Solution of the Dispersion Relation 

The reconstruction theorem has allowed us to decompose the form factors into functions of one variable, (67). 
The nine functions of one variable are given unambiguously by the Omnes solutions (85). The hat functions 
appearing in the dispersive integrals are given by angular integrals of the nine functions of one variable and 
link these functions together. Therefore, we face a set of coupled integral equations, parametrised by the nine 
subtraction constants a^®, b^°, ... as defined in (85). In this section we discuss a method for solving these 
equations numerically. We will assume an asymptotic behaviour of the form factors that requires only n = 2 
subtractions. 

A crucial property of this set of equations is that they are linear in the subtraction constants. Any solu¬ 
tion can be written as a linear combination of nine basis solutions. Our main task is therefore to determine 
numerically these nine basis solutions. 

So far, the invariant squared energy of the dilepton system, se, has been treated as an external parameter. 
On the one hand, it appears in the definition of the hat functions. On the other hand, the subtraction constants 
have an implicit dependence on si. To make this dependence explicit we write the form factors as: 

9 

X{s,t,u) ='^a,{se)X,{s,t,u), (86) 

where X € {F, G}, {ai}i = {a^° ,b^°,... ,b^°} and where Xi denotes the basis solution with Uk = Sik, 
k G {!,..., 9}. If Se is allowed to vary, the ‘functions of one variable’ become actually functions of two variables, 
Mo{s,se), ... 

Our strategy is as follows. We determine the basis solutions by a numerical iteration of the coupled integral 
equations. Each basis solution is a function of s, t and u, where s-|-f-|-u = So= -|- 2M^ + se, or 
equivalently a function of s, se and cos 6. Since se is a fixed external parameter in the integral equations, the 
iterative solution has to be performed for each value of se separately. Once the basis solutions are computed, the 
subtraction constants (or rather functions) have to be determined by suitable means, such as a fit to data, the 
soft-pion theorem and xPT input. As the dependence on se has been found to be rather weak in experiments, 
the subtraction functions can be well approximated by a low-order polynomial in se- 

In summary, we need the nine basis solutions for a set of values of se, so as to allow us to calculate them 
for any value of se by interpolation. Again, since the dependence on se appears to be rather weak, we will need 
only a low number of values of se- 

To calculate each of the basis solutions we use the following iterative procedure: 

1. set the initial value of the functions Mg, ... to Omnes function x subtraction polynomial (the polynomial 
is in fact either 0 or a simple monomial with coefficient 1 for a particular basis solution); 

2. calculate the hat functions Mg, ... by means of angular integrals of the functions Mg, ... ; 

3. calculate the new values of the functions Mg, ... as Omnes function x (polynomial + dispersive part), 
where in the dispersion integral the hat function calculated in step 2 appears; 

4. go to step 2 and iterate this procedure until convergence. 

It turns out that this iteration converges quickly. After five or six iterations, the relative changes are of order 
10 -®. 
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3.2 Phase Input 

3.2.1 TTTT Phase Shifts 

For the pion scattering phase shifts, we use the parametrisation of [23, 24]. The solution depends on 28 para¬ 
meters that can vary within a certain range. The curve labeled as Solution 1 in figure 4 shows the central solution 
for the phase shifts as well as the error band due to uncertainty in the parameters (summed in quadrature). 



0 0.2 0.4 0.6 0.8 1 1.2 1.4 1.6 1.8 2 0 0.2 0.4 0.6 0.8 1 1.2 1.4 1.6 1.8 2 

y/s/GeV y/s/GeV 


Figure 4: tttt phase shift inputs 

Two aspects deserve special attention. First, the phase for Solution 1 is just taken constant above ^/s « 
1.5 GeV. Our derivation of the dispersion relation, however, relies on the assumption that Sq(s) —> tt, t tt 

for s —> oo. We should therefore change the high-energy behaviour of the phases such that they reach tt at 
s = A^. The exact way how this is achieved should not have an influence on the result at low energies, especially 
in the physical region of the decay. We choose to interpolate smoothly between the value of Solution 1 and tt: 

'Jo(s)sol.2 := (1 - /int(si,S2,s))(jo(s)sol.l +/int (si, S2, s)7r, 

i5i(s)so 1.2 := (1 - /int(si,S 2 ,s))(j|(s)sol.l +/int (si, S 2 , s)7r, 

f 0 if s < Si, (87) 

/int(si, S 2 , s) := < if Si < S < S 2 , 

1 if S 2 < s. 

Figure 4 shows Solution 2 with si = 68M^ and S 2 = 148M^. These values can be varied and should not have 
an important influence. 

The second subtlety is the problem of the behaviour around the KK threshold [25]: are the K (4 partial 
waves expected to have a peak or a dip in the vicinity of the KK threshold, i.e. do they rather behave like the 
strange or the non-strange scalar form factor of the pion? The answer to this question could be obtained from 
a coupled-channel analysis of the K £4 amplitude, which however goes beyond the scope of this paper. In case 
of a dip we have to modify the phase such that it follows Sq(s) — tt above the KK threshold. The third solution 
shown in figure 4 is given by 

150 ( 5 ) 301.3 := (1 - /int(si,S 2 ,s)) ((5o(s)sol.l - /int(si,S2,s)7r) -f/i„t(si, S 2 , s)7r, (88) 

with Si = and S 2 = si -I- 

The Solution 4 in figure 4 corresponds to Solution 2 but with si = 4M|- and S 2 = si -I- M^. 

As the question of the correct behaviour around the KK threshold is not easy to answer, we declare Solution 3 
as the ‘central’ one and use all the other solutions to determine the systematic uncertainty. 
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Figure 5: K-k phase shift inputs, isospin I = ll2 



s/i/GeY Vt/GeV 

Figure 6: K-k phase shift inputs, isospin I = 3/2 


3.2.2 Ktt Phase Shifts 

For the crossed channels, we also need the Ktt phase shifts as an input. We use the phase shifts and uncertainties 
of [26, 27], but add by hand a more conservative uncertainty that reaches 20° at t = 15QM^. For the very small 
phase we just assume a 100% uncertainty. These phase solutions are shown in figures 5 and 6 as ‘Solution 1’. 

In the derivation of the dispersion relation, we assume that the Ktt phases go to zero at high energies. 
We implement this by interpolating smoothly between Solution 1 and zero with These modified 

phase shifts with ti = 150M^ and t 2 = 250M^ are displayed as ‘Solution 2’ in figures 5 and 6. The difference 
between ‘Solution 1’ and ‘Solution 2’ is taken as a measure of the systematic uncertainty due to the high-energy 
behaviour of the Ktt phases. 


3.3 Omnes Functions 


In a first step, the six Omnes functions are computed, defined by 


Of (s) := exp f - r , 

K'^Jso [s'-s-ie)s' ) 


where sq denotes the respective threshold. We show only the results for the tttt Omnes functions, see figures 7 and 
8 . In the case of Og, the Omnes function is computed for the phase Solution 3 - the corresponding uncertainty 
is obtained by summing in quadrature the variations generated by the uncertainties of all 28 parameters. 
The differences, appropriately weighted, are summed up in quadrature to give the error band. For the phase 
Solutions 1, 2 and 4, only the central curve is shown. Note that the differences between the various high-energy 
phase solutions are much larger than the error band due to the phase parameters. However, at low energy these 
differences are well described by polynomials and can be absorbed at low energies by the subtraction constants 
of the dispersion relation. This implies that the uncertainty generated by the unknown high-energy behaviour 
of the phase shifts will be moderate. 
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Figure 7: tttt 5-wave Omnes function, isospin 1 = 0 
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Figure 8: tttt P-wave Omnes function, isospin 1 = 1 

3.4 Hat Functions and Angular Projection 

During the iterative solution of the dispersion relation, the hat functions have to be computed by means of 
angular averages. Since the hat functions appear in the integrand of the dispersive integrals, they have to be 
known just on the real axis above the threshold of the respective channel. 

In the Si = 0 case, the calculation of the angular integrals is straightforward. The functions Mq, ... need 
to be computed on the real axis, also for negative values of their argument. As described in section 2.4.5, a 
subtlety arises in the case yf 0: in the calculation of the s-channel hat functions, we have to know the angular 
integrals of the t- and rt-channel functions Nq, .... In the region {Mk — < s < {Mk + the angular 

integration path extends into the complex t- or u-plane. Therefore, the t- and u-channel functions Ag, ... have 
to be computed not only on the real axis but also in the complex plane. Since the region where this happens is 
much below the t- or ■u-channel cut, we have two options how to perform this: 

• integrate on a straight line in the complex t- or u-plane. The functions Ngit), ... have to be known in an 

egg-shaped region of s^-dependent size. The egg lies within — Mk^/h < Re(t) < + Mk^/h- The 

functions Ao(t), ... can be computed on a two-dimensional grid covering this egg and then e.g. interpolated 
with a 2D spline. 

• since the functions Ao(t), ... are analytic in the region of the egg, the angular integration path can be 
deformed to lie always on the border of the egg. Therefore, the functions Ao(t), ... only have to be 
computed on points lying on this border (in addition to the real axis) and ID interpolation methods can 
be applied. 

The first method is more straightforward, the second needs less computing time. The second one requires a 
change of variable that we briefly describe. 
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Figure 9: Angular integration contours for se = 


We want to compute the angular integral 

(z”X)t,(s) = ^J dzz^X{t{s,z)), 

where e.g. X = Nq and 

t{s,z) = ^ (^0 “ ^ “ ^Ke{s)<^^{s)z^ . 

The square root of the Kallen function is defined by 


+ Wki{s)\ 


s < {Mk - v^)^, 


47(s) = <1 +*|A^/"(s)| (Mk - < s < (Mk + 7^)2, 

-|Ak/"(s)| (Mk + v^)2<s 

and the critical region is s_ < s < S-|-, where we define 

s± := (Mk ± vT^)^- 


(90) 


(91) 


(92) 


(93) 


In this region, the angular integration path in the complex t-plane runs from t- := t(s, —1) to := t(s, 1). 
Due to the analyticity of the function X{t), the straight contour can be deformed along the border of the egg, 
either to pass ti := t{s-,z) or t 2 ■= t{s+,z), see the two plots in figure 9. Defining 


Zsit) = 


1 




(So - s - 2t), 


(94) 


we rewrite the angular integral as a complex integral: 


1 


dzs 


= o -^dtz:{t)X{t) 

^ J 


dt 

1 _ [' 

1 / 


dtz:{t)x{t) 


(95) 


A«(s)ct^(s) 


or equivalently 


= - 


Ak^^(s)ct^(s) yt- 


dtz^mit)- J\tz^it) X it)] 


dtz:it)xit)- J\tz^it)x{t)\ 


We parametrise the border of the egg by the following curves: 

tj(e) := t(7±l) = ^ (So ? e 


(96) 


(97) 
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zntti0)xitt{0) , 
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^«(s)0'7r(s) \ds 
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(0 1 / ^ d{\^ll{i)aA0) 


di 
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1 / ^ _ 2M4 m2 - Ml {AM^se + e) + (s, - e) {2MU^ + s,) - 
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= 2 '-'^ 




ZsitfiO) = 1/2 -— • 


Note that 


zMiO) = -zsit-{or, 
tt{0 = t7{0l 
dttiO fdtjiO' 


d^ 


d^ 


and hence, due to the Schwarz reflection principle 

xitt{o) = xit:{or. 

Therefore, the function X has to be computed only on the ‘upper half-egg’: 


{z^X)t^ = 


1 


(^)^7r(^) ds — 


d^ 


dts iO 
d^ 


z:it-{o)x{t:io) 


-{-ly 


dts (0 

d^ 


z:{t:io)xit:{o) 


(98) 


(99) 


( 100 ) 


( 101 ) 


( 102 ) 


(103) 


or 


{z^xy^ 


1 

(s)0-7r(s) 



(^d^z7it-iO)xit-iO) 

- (-1)” (^^^r(t;(e))^(i;(0)) ^ • 


(104) 


Although both descriptions are valid in the range s_ < s < s+, one may choose to use the first in the region 
S- < s < Sm and the second in the region Sm < s < Sj^, where Sm lies somewhere in the middle of s_ and s+, 
e.g. Sm = (s- + s+)/2. The motivation to do so is to avoid numerical instabilities: the integral from s_ to s 
with s —)■ s+ must tend to zero to give a finite value for the hat function. The integral over the whole half-egg, 
however, accumulates a numerical uncertainty. 


30 



3.5 Results for the Basis Solutions 


We have now all the ingredients to compute the nine basis solutions of the dispersion relation. The final result 
will be a linear combination thereof. In section 4, we will describe how to determine this linear combination. 
We will fit experimental data on the partial waves defined by 

Fs{s,Si) = (^Mois,Si) + Mo{s,st)^ ^ 

Fp{s,Si) = (^Mi{s,Si) + (105) 

Gp{s,st) = + 

The figures 10 and 11 show the partial waves of the basis solutions in the case si = 0. They are computed with 
the phase solutions that reach the asymptotic values of tt in the case of the tttt phases and 0 in the case of the 
Ktt phases. For (5 q, the solution with the drop around the KK threshold is used. The figures illustrate what 
can be learnt also from the definitions (85) and (105): the data on the partial wave Fs will constrain mainly the 
subtraction constants appearing in Mq, the data on Fp mainly the constants in Mi and the data on Gp mainly 
the constants in Mi. An exception is the constant through the hat functions, it is constrained by the data 
on all partial waves. 



Figure 10: S-wave of the form factor F for the different basis solutions for si = 0 




Figure 11: P-waves of the form factors F and G for the different basis solutions for sr = 0 

Besides experimental data on the partial waves, we will also use two soft-pion theorems as an additional 
source of information. Table 1 shows the values of {F — G)(M,^, M|., M^) and {F + G){M!^, M^, M'^) for the 
basis solutions. The first soft-pion theorem, which implies {F — G){M‘^,M‘^,M‘^) k, 0, constrains mainly a 
linear combination of and 6^®. 
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basis solution 

(P — G)sppi 

[F + G)spP2 

qMo 

1.06 

1.05 

bMo 

0.08 

0.09 

^Mo 

0.03 

-0.01 


-1.03 

0.93 

bM^ 

0.05 

0.11 


-1.07 

1.02 

bM. 

-0.05 

0.09 

cMi 

-0.10 

-0.01 

bNo 

1.62 

-0.01 


Table 1: Values of the two relevant combinations of the form factors F and G at the soft-pion points, computed for the 
basis solutions. 


4 Determination of the Subtraction Constants 

In the previous chapter, we have described how to solve numerically the Omnes dispersion relation for the form 
factors F and G. The solution is parametrised in terms of the subtraction constants , .... The next task 
is now to determine these subtraction constants in order to fix the parametric freedom. We use three different 
sources of information for the determination of the subtraction constants: 

• experimental data on the form factors, 

• the soft-pion theorem, providing relations between F, G and the vector form factor, 

• input from yPT. 

The soft-pion theorem (SPT) is valid up to corrections of 0{M^) and hence can be considered as a strong 
constraint. From the two high-statistics experiments NA48/2 and E865 we have data on the S- and P-waves 
of the form factors. Although these experiments have achieved impressive results, the data alone does not 
determine all the subtraction constants with satisfactory precision. Therefore, we use chiral input to fix some 
of the subtraction constants that are not well determined by the fit to data. 

In the following, we describe in more detail what data we use for the fits and how these fits are performed. 
We were provided with additional unpublished data from the E865 experiment and include the data sets of 
NA48/2 that became available only recently as an addendum to the original publication [9]. Therefore, our fits 
include the maximal amount of experimental information on the K 14 form factors F and G that is currently 
available. 


4.1 Experimental Data 

The NA48/2 experiment defines the partial wave expansion of the form factors as 


F = FsF^^ -f FpF^^ cose + ..., 
G = GpF^’’ + ... 


and further defines the linear combination 

Gp := Gp + 

For us, it is convenient to define the partial wave 


^ a^PL 


Ml 

F--fW F 


(106) 


(107) 


(108) 


In our former treatment of the form factor Pi [12, 13, 14], it was most convenient to use the data on Fg and 
Gp (which corresponds to the P-wave of Pi). Now that we describe both form factors P and G, we prefer to fit 
the three partial waves Pg, Fp and Gp. 


32 



The comparison with our definition of the s-channel partial-wave expansions 


allows us to identify 



(109) 


( 110 ) 


The phase shifts are just given by the tttt phases that we use as input. With (80), we find the fitting equations: 

Fsis,Si) = (^Mo(s,S() + Mo{s,Si)^ ^ 

Fpis,st)= (Mi(s,Sf) + Mi(s,s^))e-*'^^(^), (111) 

Gp{s,si) = (Mi(s,s^) + Mi(s,s,)) 


The NA48/2 collaboration has performed phenomenological fits of the form [6, 9] 


Fs{s,se) 2 ,11 4 , f'e se 

fs 


Fp{s,se) 

fs 


fp 

fs' 


Gp{s, sg) 9p , 9p 2 

fs fs fs^ ' 


( 112 ) 


where — 1. In a first step, only the normalised coefficients were measured [6]. In a second step, the 

normalisation fs was determined from the branching ratio measurement and a phase-space integration, using 
the parametrisation (112) and the fitted normalised coefficients [9]. 

However, one should note that from (110) it follows that Fp has to vanish at the tttt threshold like ~ \fq^- 
The phenomenological fit (112) of [6, 9], which assumes Fp to be constant in gives a wrong threshold 
behaviour. We have not tried to estimate its influence on the determination of the normalisation fg. For our 
purpose, we find it convenient to work with Fp, which does not contain kinematic prefactors. 

Because all the basis solutions use the same tttt phase as input, the real quantities Fg, Fp and Gp are still 
linear combinations of the corresponding quantities computed with the basis solutions. Note that the partial 
waves can be negative, i.e. one really has to rotate the tttt phase away and not just take the absolute value. 

For our fits, we use the experimental values of NA48/2 [6, 9] and E865 [7, 8] on the partial waves. Some 
remarks on these numbers are appropriate. 

• Originally, the published NA48/2 data consisted of 10 bins in s-direction. Very recently, a two-dimensional 
data set on Fg{s, Si) has become available (addendum to [9]): in this set, not only a single bin but up to 
10 bins are used in s^-direction. 

• The barycentre values of Sf, for the original 10 bins of NA48/2 also became available in the addendum to 
[9]. A value of S( could also be extracted from the relation (107) between Fp, Gp and Gp [28]. However, 
this value does not agree with the barycentre. 

• We compute the value of Fp with (108) using the values of Fp and the barycentre values of s and se. 

• There is a discrepancy between [6] and [9]. The statistical and systematic uncertainties for Fg in the 
NA48/2 data have to be calculated from the normalised coefficients in [6]. The correct uncertainties are 
also listed in the addendum to [9]. 
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• The published values of Fg in the 10 bins of NA48/2 have been normalised in such a way that a fit of the 
form (112) with /e = 0 results in Fs{0,0)/fs = 1, although a non-zero value of /' has been obtained from 
a fit to the two-dimensional data set. In order to take the s^-dependence consistently into account, we 
have to increase the values of Fg by 0.77%. 

• The E865 experiment has assumed in the analysis that the form factors do not depend on Si. The values 
of Si for each bin were not published.^ 

• The E865 experiment only provides data on the first partial waves Fg and Gp. 

• The E865 papers [7, 8] include the fully correlated error of the normalisation of 1.2% in their systematic 
errors (added in quadrature).^ It needs a special treatment for unbiased fitting. 

In the data analysis of both experiments, radiative corrections have been applied to some extent. More 
reliable radiative corrections based on a fixed-order calculation [28] can be applied a posteriori at least to the 
NA48/2 data. Furthermore, neither the E865 nor the NA48/2 experiment has corrected the isospin-breaking 
effects due to the quark and meson mass differences. The calculation of [28] also allows for their correction. The 
resulting numbers are given in appendix D. We add the uncertainties of the isospin corrections (without the 
higher order estimate) in quadrature to the systematic errors. The one-dimensional NA48/2 values also include 
the mentioned correction of the normalisation of Fg by 0.77% due to the s^-dependence. 

In addition to the statistical and systematic errors, we take into account the correlations between Fp and 
Gp of the NA48/2 data, which also became available with the addendum to [9[. There are however several 
correlations that we neglect, either because they are not available or because we assume them to play a minor 
role. These include the bin-to-bin correlations of the P-waves and the correlations with the 5'-wave. We also 
neglect the correlation due to the isospin-breaking corrections and correlations between the two experiments 
due to external input. We do not expect any of these neglected correlations to significantly affect our fits, but 
of course it would be better to check. If the complete set of experimental correlations will become available, it 
will be possible to do that. 


4.2 Soft-Pion Theorem 

In addition to the experimental input on the partial waves, we use the soft-pion theorem (SPT) [29, 30] as a 
second source of information to determine the subtraction constants. 

There are two different soft-pion theorems for Ki 4 , depending on which pion is taken to be soft. If the 
momentum pi of the positively charged pion is sent to zero, the Mandelstam variables become s = M^, t = M'^, 
u = Si- Since the SPT is valid only at 0{M^), we set u = Si+M^, such that the relation s+t+u = M'^+2M^+Si 
remains valid and one does not need to worry about defining an off-shell form factor. 

The first SPT states [12]: 

FiMl Mi, m 2 + s,) - G(Mi Mi, M^ + s,) =0{M^). (113) 


If the momentum p 2 of the negatively charged pion is sent to zero, the Mandelstam variables become s = M^, 
t = Si, u = Mi. We set t = Si + Mi 

The second SPT gives a relation to the Kis vector form factor: 


F{MlMl + Si, Mi) + G(Mi Ml 


se, Mi) 


y/^Mx 

F^ 


U{Ml + si)=0{Ml). 


(114) 


At leading order in yPT, the SPTs are fulfilled exactly, i.e. the right-hand sides of the equations (113) and 
(114) vanish, at NLO and NNLO, there appear 0{Ml) corrections. 

Numerically, it turns out that the first SPT is fulfilled to a higher precision than the second SPT. At NLO, 
the correction to the first SPT is about 0.4% for si = 0, the second SPT gets a correction of 2.0% if f+{Ml) 
is used. If we make the arbitrary replacement /^.{Ml) i—>■ /+(0), again an 0{Ml) effect, the deviation in the 
second SPT increases to 4.9%. This confirms that the size of the observed deviations from the SPT is natural. 

At NNLO, the corrections become slightly larger.^ If the 0{p^) LECs Gf are all put to zero and = 0 as 
well, the first SPT is fulfilled at 1.0%, the second at 4.4% with /_|_(M2) or 7.6% with /+(0). If the Gf parts are 

^We thank Peter Truol and Andries van der Schaaf, who performed a new analysis of the Brookhaven data in order to extract 
the barycentre values of . 

^We thank Stefan Pislak and Peter Truol for this additional unpublished information. 

^We thank Johan Bijnens and Ilaria Jemos for providing the CH—h implementation of the NNLO expressions. 
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replaced by the estimates of [11, 31] (resonances estimates in the case of K 14 ), the accuracy of the first SPT is 
1.5%, the one of the second SPT 5.4% using /+(M^) or again 7.6% using /+(0). 

We use the size of the NNLO corrections to the SPT as an estimate of the tolerance that we allow in the 
fits when using the SPTs as constraints. 


4.3 Fitting Method 

In the following, we describe how we perform the fit. Basically, we have to deal with a simple linear fit. The 
only subtlety is the fact that the data contains a fully correlated uncertainty of the normalisation, which is a 
multiplicative quantity. The fact that we use two experiments with different normalisation errors asks for a 
special fitting method to avoid a bias [32, 33[. We apply the ‘to-method’ of [33[. 

First, we construct a covariance matrix for the observations as follows. 

• For all the partial-wave data that we want to fit we construct the covariance matrix with the squared sta¬ 
tistical errors on the diagonal and the statistical covariance between the P-waves as off-diagonal elements. 

• We add the uncorrelated systematic errors, which do not contain the error of the normalisation, in quadra¬ 
ture to the diagonal entries. 

• We may or may not include the two soft-pion theorems as additional observations. If we do so, we take 
e.g. F—G at the first soft-pion point (SPP) and F+G at the second SPP as observations. As uncertainties, 
we take a value typical for the deviation in yPT at NNLO, e.g. 1% or 2% of the LO value of F for the 
first SPT and a few percent of •\/2Mif/Pn./+(0) for the second SPT. 

• We add the errors of the normalisation to the covariance matrix, which are in block-diagonal form for the 
data of the two experiments: 

(cov)y = (rel.cov.)ij -|- (norm.cov.)*^, (norm.cov.)^^ = Aj f{s\ (115) 

where A/ denotes the error of the normalisation for experiment I. A is the index of the experiment (1 or 
2) corresponding to the data point i and /(s^s^) is the value of the fitted partial wave. In a first step, 
this value has to be computed under the assumption of some starting values for the fit parameters. 

The fit requires then an iteration. One has to minimise the error function defined by 

X" = v'^Pv, (116) 

where v is the vector of the residues, i.e. the differences between the observations and computed values. P is 
the inverse of the covariance matrix constructed above: P = (cov)“^. The minimum of the function can be 
either found with some minimisation routine or, since the fit is linear, directly with the explicit solution 

par = {A^PA)-^A^PO, (117) 


where O is the vector of observations and 


df{s\s\) 

9param^- 


(118) 


is the design matrix to be determined with the values of the basis solutions. 

With these new values for the fit parameters, one again computes the new covariance matrix (the contribution 
for the normalisation changes) and iterates this procedure. It turns out that only very few iterations are needed 
to reach convergence. 

If we do not want to determine a parameter through the fit but fix it beforehand to a non-zero value, we 
have to subtract the fixed contribution from the observations O, such that O is purely linear in the parameters 
and contains no constant contributions. 

In the above discussion, we have not specified what we use as fit parameters. One option is to fit the 
subtraction constants. Since we want to include an s^-dependence in the subtraction constants, we write e.g. 




(119) 


where a^“, a^“, ... are now the parameters collected in the above vector ‘par’. Another option is to use the 
matching equations to yPT, which provide a linear relation between the subtraction constants and the LECs 
we are interested in, and perform the fit directly with the LECs. 
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4.4 Matching to xPT 

The final goal of this treatment is the determination of low-energy constants of xPT. Instead of fitting directly 
the Kii data with the chiral expressions, we use the dispersive representation as an intermediate step. The 
dispersion relation provides a model-independent resummation of final-state rescattering effects. Therefore, 
we expect that even the most important effects beyond 0{p^) are included in the dispersion relation. Of 
course, in order to extract values for the LECs, one has to perform a matching of the dispersive and the chiral 
representations. This can be done e.g. on the level of the form factors [12, 13, 14]. Since the dispersion relation 
describes the energy dependence, the matching point can be outside the physical region, i.e. even at lower 
energies, where yPT is expected to converge better. 

Here, we follow an improved strategy for the matching: we match the dispersive and the chiral representations 
not on the level of form factors but directly on the level of subtraction constants. Since the decomposition (67) 
is valid up to terms of 0{p^), the one-loop and even the two-loop result can be written in this form, which 
allows us to extract a chiral representation of the subtraction constants. This procedure has the advantage that 
the matching is performed for each function of one variable Mq{s), ... at its subtraction point, i.e. at s = 0, 
t = 0 and u = 0, where indeed the chiral representation is expected to be reliable. 


4.4.1 Matching Equations at 0{p'^) 

4.4.1.1 Reconstruction of the yPT Form Factors 

Let us start by reconstructing the NLO form factors in the standard dispersive form (68). 
The LO xPT form factors are given by 


T^lo = Glo 


Mk 

V2F^' 


With the partial wave projections (58), we find 


/i"°(s) 




Mk 

Mk M^PL 
V2F^ 2X2 
Mk 
V2F^’ 


The isospin 1/2 form factors (20) are given by 


p(i/2) _ Mk Mi/2) _ V^Mk 

“ V2F^" “ F^ 


Hence, the partial waves in the crossed channels (60) are 


f(l/2) 

JO,LO 


it) 


f(l/2) 

Jl,LO 


it) 



Mk 5Aktt — 5 t 

V2F. 4t ’ 

Mk 3Mf,{M^-s,-t) 

y/2F^ 

3Mk 


.(3/2) 

/0,LO 


(u) 


Mk 

V2F,’ 


.(3/2) 

/l,LO 


(u) = 0, 


ffpLO = 0- 


The TTTT-scattering amplitude can be written as [2] 


u) = 34(s, t, u) + A(t, u, s) -I- A{u, s,t), 
r^^^(s,f,u) = A{t,u,s) — A{u,s,t), 


( 120 ) 


( 121 ) 


( 122 ) 


(123) 


(124) 
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where at LO 


t, u) = 


s-Ml 


FS ' 

TT 

The Mandelstam variables for tttt scattering satisfy 

s + t + u = 4M^, 

t = -2q^{l - z), 

where z = cosd. Hence, the tttt partial waves are 


(125) 


(126) 


<lo(5) = I 

4,lo(s) = ^ J ^dz zTIII (s, z) = 


s — 4M^ 


(127) 


The ifTT-scattering amplitude is given by [34] 


and at LO 




r(3/2)(s,t,w) = —(Mi + M^-s). 


(128) 


(129) 


Of course, the Mandelstam variables satisfy here s + t + u = 2M^ + 2M^. The partial waves are given by 


^ofLo(s) = ^ - MM], + Ml) - 5AM) , 


^i)lo('®) “ SsF'^ ~ + Ml) + 3A^„.) , 


^o!lo('®) — 2 ^Mk + Ml - s), 

+3/2 


(130) 


tl'loM = o. 

Using the unitarity relation for the partial waves, we can now easily construct their imaginary parts at 
NLO: 


ImfM^is) 

ImgM^is) 

I™//,NLO (^) 

(^) 


1 1 
2/+ 13^ 
1 1 


21 + 1 327r 

^ ^ ^kM ) +/♦ 

2 / + 1 IGtt t 

^ ^ ^Ktt (^) +/♦ 

2 / + 1 IGtt t 


<^^is)tiXois)fMis), 


(131) 
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hence 


Im/o,NLo(0 

W1nlo(0 

Im5i^NLo(0 

Wd^NLo(“) 

^^fi%Loiu) 


1 ,,Mk{2s-M^) 

V2F3 ’ 

1^ Mk{s- 4 .M^)MIPL 
3 y2i^3 2X2 > 

1 1 Mk{s-AMI) 

-a^(s)- 


3 327r y2i^3 

1 Mk 

167r t 32V2PF3 

1 Ad/^(i) Mk 


(5^2 - 2 t(M|. + m 2 ) - 3A^,) (3A;f, - 5t), 


IhTT t SV2tF^ 

1 Mk 

IhTT t 16v^iF3 

1 Add^(z.) Mk 

167 r u 2y/2F^ 
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0 . 


(3i2-6i(M|.+M2) + 3A2^J 


M4 (M2 _ 




(3t2-6i(M|. + M2)+3A^J, 


(M^ + m2 - m), 


( 132 ) 


By inserting these imaginary parts into the dispersion integrals in ( 68 ), we can reconstruct the NLO form 
factors. For the comparison with the explicit loop calculation, we rewrite the dispersive integrals in terms of 
loop functions (see appendix A): 


Mq^^{s) — Wq nLO + + 2y2F3 ~ (B-KTris) — i?7r7i-(0)) + M^s , 

Mi^^O(s) = m?^NLO> 

MJ^^°(s) = nLO + ^i.NLO)^^ + Q^p3 ~ 4M2) (i?,r7i-(s) — .B7r7i-(0)) + dM^S , 

+ 32^ ( (-25^ + 5(5M|. - M^)) (B;,^(t) - BkAO)) 

+ (i(3A^i - 7M2) - 3Al^) (^BK^t) - BkAO) -tB'^AO) - | ^".(0)) 

- 5t(5M| - Ml)B'j,^{Q) + , (I33) 

^NLO(^) ^ 

jYNLO(^) _ ___^_ ^^2 _ 2t{M‘^ + Ml) + A|-,r) {BKirit) — Bktt{^)) 

+ i (2t(Mi + m2) - A^J ,B),,(0) - ^ B^.(O)), 

= 2^ ( i^K + Ml-u) {BkAu) - BkA^)) - {Ml + Ml)uBl^{Q)) , 

R^^^{u) = 0, 

R^^^{u) = 0. 


We can now compare this expression with the one-loop calculation [35, 36, 10]. As in our dispersive treatment, 
we only consider tttt intermediate states in the s-channel and Ktt intermediate states in the crossed channels, 
the KK and 7777 loops in the s-channel and the Krj loops in the f-channel have to be expanded in a Taylor series 
and absorbed by the subtraction polynomial. The comparison of the dispersive representation with the loop 
calculation then allows the extraction of the 0{p‘^) values for the subtraction constants. 
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Note that the only contributions that we neglect when writing the loop calculation in the dispersive 

form are the second and higher order Taylor coefficients of the expanded loop functions of higher intermediate 
states {KK, 7777 and Krj). The result for the 0{p^) subtraction constants can be found in appendix E.l. 


4.4.1.2 xPT Form Factors in the Omnes Representation 

The reason why we do not use the standard dispersive form ( 68 ) for the numerical solution of the dispersion 
relation but rather the Omnes representation (85) is mainly the separation of final-state rescattering effects: 
the Omnes function resums the most important rescattering effects. The remaining dispersive integrals take the 
interplay of the different channels into account. 

It is therefore desirable to perform the matching to yPT not on the level of the standard dispersive form 
but directly with the Omnes representation. This should avoid mixing the final-state resummation with the 
determination of the LECs. 

However, it is not possible to write directly the yPT representation in the Omnes form, because the chiral 
expansion of the phase shifts does not have the correct asymptotic behaviour. At LO, the phases grow linearly, 
hence the Omnes dispersion integral (83) is logarithmically divergent. Therefore, we subtract the dispersion 
integral once more: 


2 /*oo 


= exp [ - 

TT 


> So 
/*oo 




«(»') 


TT (s' - s - ze)s'^ 


,ds' 


(134) 


<5(s') 
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TT Jao {s' - S- 


Tds' I . 


oj is divergent if evaluated in yPT. Let us postpone the determination of this constant for a moment. 

Let us now use the Omnes representation to reconstruct the NLO result for the form factors. At LO, the 
functions of one variable are simply given by 


mLO(s) = MLO(s) = ^, 

ML°(s) = N^'^{t) = = iVi^°(f) = i?^°(w) = i?i°(u) = Ri°(u) = 0. 

We start by calculating the hat functions at LO: 


MLO(s) 

fV“(<) 

i?LO(u) 
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Mk iMl 
2t ’ 
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Further, we need the phase shifts at LO: 
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(137) 
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We expand the Omnes representation (85) at NLO: 
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If we further expand these expressions chirally and neglect higher orders, we obtain (note that only a^® and 
a^i do not vanish at LO): 
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Next, we insert the LO phases and hat functions: 
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We see that the form of the Omnes representation is completely equivalent to the standard representation, apart 
from the presence of the additional subtraction constants c^°, and which also need to be determined. 
We expand the t-channel Krj integrals up to linear terms in t and find: 
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(142) 


The constants Wg and wj cannot be evaluated with the chiral phases. If we evaluate them with the physical 
phases, this leads to exactly the same matching equations for the determination of the as if we would match 
the Taylor expansion of the Omnes representation with the Taylor expansion of the chiral result. Note, however, 
that the expressions obtained for c^°, b^^ and are different. E.g. for 6^0 the chiral expansion leads to 
6^0 = 0 while a Taylor expansion of the dispersion relation would require b^^ = —to® nlo^i i^)Ml, where 0} 
is the derivative of the Omnes function calculated with the physical phases. Of course, the difference is a higher 
order effect in the chiral counting. As higher order effects can be important if due to final state rescattering, we 
would not like to intermingle them with the matching of the subtraction constants. The matching on the basis 
of Taylor coefficients would require the linear term of Mi{s) to vanish exactly, while the matching based on the 
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chiral expansion of the dispersion relation gives a non-zero linear term in Mi(s) due to the Omnes function - 
this is important information which we wish to make use of in our fits. 

4.4.2 Matching Equations at 0{p^) 

4.4.2.1 Decomposition of the NNLO Form Factors 

In the following, we describe the decomposition of the two-loop result such that the matching can be performed 
at NNLO. Since the NNLO chiral result has a different asymptotic behaviour than the NLO result and our 
numerical dispersive representation, we have to use the representation (69), which uses a different gauge and 
more subtractions than (68). 

The imaginary parts of the K 14 partial waves at NNLO could again be reconstructed using the unitarity 
relations, e.g. 

Im/;^^^°(s) = 2 ; _|_ 1 (^f.Lo('S)/i"°(s) + ^^f,NLo(s)/^°(®) + • (143) 

However, instead of proceeding as for NLO, it is more straightforward to decompose the two-loop result directly 
into functions of one variable, then to impose the gauge condition and extract the Taylor coefficients of the 
functions of one variable. 

The two-loop result for the form factors F and G was computed in [11]. We have the full expressions in 
form of a CH—h program at hand."^ It has the following structure: 

t, u) = -k t, u) -k t, U) 

+ t, u) + xNnlO(s, t, u) + xNnlO(s, t, u) (144) 

+ t, u) + t, u) + t, u), 

where X G {F, G} and the different parts denote the following: 

• X^^^: NLO polynomial containing the LECs L^, 

• X^LO. NLO loops, 

• NNLO polynomial containing the LECs C^, 

• NNLO part containing L'^ x and x loop, 

• NNLO two-loop part without vertex integrals, 

• NNLO vertex integrals in the s-channel, 

• NNLO vertex integrals in the t-channel, 

• Xyp'"®: NNLO vertex integrals in the u-channel. 

In appendix E.2.1, we perform the explicit decomposition of the two-loop result into functions of one Man¬ 
delstam variable according to (67) and (69) and evaluate numerically the subtraction constants. 

4.4.2.2 NNLO Form Factors in the Omnes Representation 

As we already pointed out for the NLO matching, it is desirable to use the Omnes representation rather than 
the standard dispersion relation for the matching and the determination of the LECs. Let us therefore derive 
the matching equations at NNLO in the Omnes scheme. 

We have to use the second gauge for the decomposition of the NNLO representation (69). As a starting 
point, let us find the NLO Omnes subtraction constants in the second gauge. In the first gauge, we found 
i?NLO ^ ^NLO ^ 0 , hence 

CnLO = ^p3 {i^K + i 

Ri _ Li?i _ r, 

“NLO “ “nLO “ 

“^We thank Johan Bijnens and Ilaria Jemos for providing the CH—h implementation of the NNLO expressions. 
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The gauge-transformation (74) is then defined by 
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At NLO, the shifts in the subtraction constants (C.5) are therefore given by 
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When studying now the Omnes representation at NNLO, we notice that the asymptotic behaviour of the 
phases at NNLO is even worse than at NLO, hence we have to subtract the Omnes function three times: 
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u) and Co are both divergent if evaluated in yPT at NNLO, hence we will use the physical phases to determine 
them. 

In the case of the NLO matching, we have derived the relation between the standard and the Omnes 
subtraction constants (142) by comparing the Taylor coefficients of the chirally expanded Omnes representation 
with the Taylor coefficients of the standard dispersive representation. Although it is instructive to understand 
the chiral expansion of the Omnes representation, a shortcut can be taken. Note that the chiral expansion 
and the Taylor expansion are interchangeable. Therefore, we easily obtain the relations between the standard 
subtraction constants TOq, ... and the Omnes subtraction constants a^°, ... by chirally expanding the Taylor 
coefficients of the Omnes representation (C.l) and comparing it with the Taylor coefficients of (69). 
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This leads to the following relations between the relevant subtraction constants: 
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The NNLO chiral expansion of the full Omnes representation can be found in appendix E.2.2 and leads to the 
same result. It can be used to identify all the imaginary parts and to connect the different dispersive integrals 
with the discontinuities of the loop diagrams. 


5 Results 

In this chapter, we discuss the results for the low-energy constants that we determine by fitting the dispersive 
representation to data and matching it to xPT. In order to understand the differences between the results at 
NLO and NNLO and the source of complications that appear at NNLO, it is useful to study in a first step the 
results of direct yPT fits. We perform direct fits at NLO and NNLO and compare our results with the literature 
before using the whole machinery of the dispersive framework matched to yPT at NLO and finally at NNLO. 

5.1 Comparison of Direct xPT Fits 

The most recent fits to Ki 4 data performed in the literature are [37]. There, a global fit is performed, taking 
into account the threshold expansion parameters of the K (4 form factor measurement of NA48/2 [9]: 

F = fs + fiq^ + ..., /,= 5.705 ± 0.035, /( = 0.867 ± 0.050 

, 9 , (150) 

G = gp +g'^q^ + ..., c/p = 4.952 ± 0.086, 5 p = 0.508 ±0.122, 

where — 1. In [37], the above quantities are fitted with the form factors at cosd = 0 instead of the 

first partial wave. In addition to the K £4 form factor data, the global fit of [37] uses many other inputs, like 
data on the different decay constants and masses, tttt- and iL 7 r-scattering parameters, quark mass ratios etc. 

We compare now different strategies for direct fits with the results of [37]. We use only K 14 data for our fits 
and therefore are only sensitive to the LECs L\, ±2 and L 3 [12]. The other LECs are taken as a fixed input. 

5.1.1 Direct Fits at 0{p'^) 

5.1.1.1 Fits of Threshold Parameters 

In order to make the connection to [37], we first perform a direct NLO fit to the NA48/2 threshold parameters 
in (150). Using cos 6 * = 0, i.e. the first Taylor coefficient of an expansion in 2 ; = cos 6 , and the LEG inputs 
±4 = 0 and the fitted value for ±5 of [37], we reproduce almost exactly the result of [37] for LJ, ±2 and ± 3 , see 
the second and third column in table 2. If we use instead the partial-wave projection (55), the fit results for 
and ±2 change a bit, as shown in the fourth column of table 2. The last column uses lattice results [38, 39] for 
the input LECs. 
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Ref. [37] 

Taylor 

PWE 

PWE 

10 ^ • LI 

0.98(09) 

0.99(09) 

1.15(09) 

1.17(09) 

lO^ • 

1.56(09) 

1.57(09) 

1.48(08) 

1.50(08) 

10 ^ • 

-3.82(30) 

-3.83(30) 

-3.82(30) 

-3.87(30) 

103 • LI 

= 0 

= 0 

= 0 

= 0.04 

lO^ • LI 

1.23(06) 

= 1.23 

= 1.23 

= 0.84 


16 

0.3 

0.3 

0.3 

dof 

5 

1 

1 

1 

X^/dof 

3.2 

0.3 

0.3 

0.3 


Table 2: Comparison of direct NLO fits to the NA48/2 threshold parameters [9]. The renormalisation scale is /r = 
770 MeV. The last column uses the lattice determination of [38, 39] for the input LECs. The uncertainties are purely 
statistical. 


5.1.1.2 Fits of the Complete Form Factor Data 

In a next step, we no longer fit the threshold expansion parameters (150) of the form factors, but the form 
factor data of NA48/2 [6, 9] and E865 [7, 8], discussed in section 4.1. The second column of table 3 shows the 
result of the NLO fit to the one-dimensional NA48/2 data without isospin corrections (but with the corrected 
normalisation of Fs to account for the s^-dependence). In the third column, isospin corrections are applied to 
the fitted data (table 10). The fourth and fifth column show the results of a combined fit to NA48/2 and E865 
data (table 11). The smaller value in the fits to the data with isospin-breaking corrections is due to the fact 
that the isospin corrections introduce an additional uncertainty in the data. 



NA48/2 

NA48/2,isd 

NA48/2 & E865 

NA48/2 & E865,is(5' 

10 ^ • LI 

0.69(03) 

0.71(04) 

0.62(03) 

0.64(04) 

10 ^ • LI 

1.88(07) 

1.80(08) 

1.79(06) 

1.70(06) 

10 ^ • LI 

-3.89(13) 

-3.93(14) 

-3.62(11) 

-3.60(12) 

10 ^ • LI 

= 0.04 

= 0.04 

= 0.04 

= 0.04 

10 ^ • LI 

= 0.84 

= 0.84 

= 0.84 

= 0.84 

10 ^ • LI 

= 5.93 

= 5.93 

= 5.93 

= 5.93 


159.4 

67.5 

199.9 

117.1 

dof 

27 

27 

39 

39 

X^/dof 

5.9 

2.5 

5.1 

3.0 


Table 3: Comparison of direct NLO fits to the NA48/2 and E865 form factor measurements. The renormalisation scale 
is /i = 770 MeV. For and L 5 , we use lattice input [38, 39], for Lg the determination of [40[. The uncertainties are 
purely statistical. 


Figure 12 shows a comparison of the NA48/2 threshold parameter fit of [37] with the result of the fit to the 
whole form factor data set (forth column of table 3). It helps to understand the difference between the fitted 
LECs in the two procedures: in the fit to the threshold parameters, the curvature of the form factor is neglected. 
Since the NLO chiral representation cannot reproduce the curvature, the data points at higher energies reduce 
the slope in a fit to the whole data set. 
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5-wave of F 



s/GeV^ 


Figure 12: Comparison of different fits for the S-wave of the form factor F: NA48/2 threshold parameter fit of [37] 
and a fit to the full data set. The (s, se) phase space is projected on the s-axis. No isospin corrections are applied. 


5.1.2 Direct Fits at 0{p^) 



Ref. [37] 

Ref. [37] 

NA48/2 

NA48/2 & E865 

NA48/2 

NA48/2 & E865 

Cl 

= 0 

BE14 

= 0 

= 0 

BE14 

BE14 

10 ^ • LI 

0.67(06) 

0.53(06) 

0.34(03) 

0.28(02) 

0.33(03) 

0.27(02) 

10 ^ • LI 

0.17(04) 

0.81(04) 

0.42(06) 

0.35(05) 

0.95(06) 

0.89(05) 

10 ^ • LI 

-1.76(21) 

-3.07(20) 

-1.54(14) 

-1.25(11) 

-3.06(14) 

-2.80(11) 

103 • LI 

0.73(10) 

= 0.3 

= 0.04 

= 0.04 

= 0.04 

= 0.04 

10 ^ • LI 

0.65(05) 

1.01(06) 

= 0.84 

= 0.84 

= 0.84 

= 0.84 

10 ^ • Lg 

0.25(09) 

0.14(05) 

= 0.07 

= 0.07 

= 0.07 

= 0.07 

10 ^ • Lg 

-0.17(06) 

-0.34(09) 

= -0.34 

= -0.34 

= -0.34 

= -0.34 

10^ • Lg 

0.22(08) 

0.47(10) 

= 0.36 

= 0.36 

= 0.36 

= 0.36 

10^ • Lg 



= 5.93 

= 5.93 

= 5.93 

= 5.93 


26 

1.0 

81.3 

128.7 

52.5 

91.2 

dof 

9 


27 

39 

27 

39 

X^/dof 

2.9 


3.0 

3.3 

1.9 

2.3 


Table 4: Direct NNLO fits for different choices of the Cl. The results of the fits of [37] are shown for comparison. The 
renormalisation scale is fi = 770 MeV. Our results are fits to the entire form factor data including isospin corrections. 
The uncertainties are purely statistical. The NLO input LECs LJ, LI, Lg and Lg are lattice determinations [38, 39], LI 
is the BE14 value [37] and Lg is taken from [40]. 


xPT at NNLO suffers from the problem that many new low-energy constants Cl appear in the 0{p^) Lagrangian. 
In Ki 4 , in total 24 linearly independent combinations of the Cl enter in the NNLO chiral representation of the 
form factors F and C. A fit of so many parameters seems out of question. We would rather like to use some 
input values for the Cl- Unfortunately, only very few of the NNLO LECs are known reliably. We could either 
use determinations of the Cl with models like the chiral quark model [41], a resonance estimate [11, 42] or the 
educated guess of [37]. These different estimates, however, do not lead to compatible results [37]. 

In table 4, we display the results of our direct xPT fits at NNLO in comparison with the results of [37]. In 
contrast to [37], we do not use the threshold parameters but the whole form factor data sets of NA48/2 and 
E865 corrected by isospin-breaking effects [28]. It turns out that even at NNLO, yPT has trouble to reproduce 
the curvature of the Fs data. We also note that the results for the fitted LECs at NNLO differ quite significantly 
from the results at NLO. 
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5.2 Matching the Dispersion Relation to xPT 

With the direct xPT fits, we have seen a number of problems: First, at NLO and even at NNLO, the energy 
dependence of the Fs form factor is not very well described. Second, at 0{p^), the appearance of quite a large 
number of additional LECs reduces the predictive power of yPT. Some input values for the Cf have to be 
assumed, as a fit of K 14 data alone cannot determine all these LECs. 

We now turn to the results using the dispersive representation as an intermediate step in the determination 
of the LECs: we fit the Ki 4 form factor data with the dispersion relation. The matching to yPT relates the 
subtraction constants of the dispersion relation to the LECs. As the dispersion relation provides a resummation 
of final-state rescattering effects, we trust that we will obtain a better description of the energy dependence of 
the form factors. However, it is clear that the matching of the dispersion relation to NNLO yPT will not be 
free of the problem related to the large number of LECs. We will alleviate the situation by including additional 
constraints on the chiral convergence in the fit. This will enable us to fit partially the contribution of the NNLO 
LECs to the subtraction constants. 

5.2.1 Matching at 0{p‘^) 

Our numerical solution of the dispersion relation (85) is parametrised by nine subtraction constants, which in 
fact are functions of se. If we use the matching at NLO to provide a chiral representation of the subtraction 
constants, we see that a^o linear in Si, while the other subtraction constants do not depend on 

S£. We therefore introduce this s^-dependence according to (119) and have to determine in total 11 parameters. 

We fit our dispersive representation to the data of both experiments, shown in appendix D. In the case of 
NA48/2, the use of the two-dimensional instead of the one-dimensional data set has basically no effect on the 
determination of the LECs L\, L 2 and Lg but gives us the option to fit the s^-dependence and therefore to 
determine also Lg. In order to test the influence of the isospin-breaking corrections, we also perform fits to data 
without isospin corrections. 

An unconstrained fit with the 11 subtraction parameters leads to a low relative of 0.77 (with 94 degrees 
of freedom, dof) for the NA48/2 data alone or 0.74 (106 dof) for the combined data set of NA48/2 and E865. 
However, the soft-pion theorems in such a fit are not well reproduced. Therefore, we chose to use the soft-pion 
theorems as constraints in the fit: the first soft-pion theorem (113) with a tolerance of 2% and the second 
soft-pion theorem (114) of 5%. These numbers are inspired by the typical NNLO deviation. In these fits the 
relative slightly increases to 0.79 (96 dof) for the NA48/2 fit and 0.77 (108 dof) for the combined fit. This 
shows that in a fit with all 11 parameters, the soft-pion theorems are not fulfilled automatically but are not a 
strong additional constraint. 

In an unconstrained fit, the result for the subtraction constants turns out to be rather unstable: the statistical 
uncertainties are large and some of the subtraction constants change drastically if the E865 data is included. 
We consider these fits of little interest and fix to an a priori value those subtraction constants that have the 
largest statistical uncertainty: these are the subtraction constants of highest order in each function and the one 
parametrising the s^-dependence in Gp, i.e. c^°, and We fix these subtraction constants to the 

NLO chiral prediction in the matching (142): while c^°, and are purely numerical, depends on 
Lg and Lg. We take those two LECs as input and iterate the fit after the matching to reach self-consistency for 
Lg (and Lg if this LEG is determined in the matching as well). 

Seven subtraction constants a^°, , b^°, b^^ and b^^ remain to be fitted to data. In the 

matching equations (i.e. (142) together with appendix E.l), the LECs L\, Lg, Lg and Lg are overdetermined. 
Hence, we have to use a second minimisation to fix these LECs. As an alternative to this two-step procedure 
(first fit to data, then matching to xPT), we can directly use the NLO chiral representation of the subtraction 
constants and perform the fit of the dispersive representation to data with the LECs as fitting parameters. As 
expected, these two strategies lead to almost identical numerical results for the LECs. 

In table 5, we show the results of the fits of the dispersion relation matched to NLO xPT. For the input 
LECs, we use lattice results [38, 39]: 


10^ • Lf = 0.04(14), 

, (151) 

10^ • L^ = 0.84(38). 

The x^ and degrees of freedom correspond to the strategy of using the LECs as fitting parameters. If we 
use the two-step fitting/matching strategy instead, the x^/dof of the fit of the subtraction constants to data 
is good: around 0.8 for the fit to NA48/2 and around 1.0 for the fit to both experiments. At the same time. 
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the relative x^/dof of the matching is bad (between 2.9 and 6.1). This is not surprising because the sum of the 
total of the two steps is approximately equal to the total in the one-step procedure, while the dof in the 
second step are drastically reduced. 

The first bracket indicates the statistical uncertainty due to the fitted data. The second bracket gives the 
systematic uncertainty. In section 5.3, we will discuss in more detail the different sources of uncertainty. 



NA48/2 

NA48/2 & E865 

NA48/2 

NA48/2 & E865 

NA48/2 

NA48/2 & E865 

Isospin corr. 

X 

X 

/ 

/ 

/ 

/ 

(TSPTI 

— 

— 

— 

— 

2 % 

2 % 

0'SPT2 

— 

— 

— 

— 

5% 

5% 

10® • LI 

0.52(02)(05) 

0.48(02)(05) 

0.54(02)(05) 

0.50(02)(05) 

0.54(02)(05) 

0.50(02)(05) 

10® • Lg 

1.00(05)(07) 

0.94(04)(07) 

0.94(05)(07) 

0.88(05)(07) 

0.94(05)(07) 

0.88(05)(07) 

10® • Lg 

-3.03(11)(07) 

-2.83(09)(07) 

-2.99(11)(07) 

-2.79(10)(07) 

-2.99(11)(07) 

-2.80(10)(07) 

10® • Lg 

4.70(40)(63) 

4.64(39)(61) 

4.51(43)(63) 

4.44(43)(61) 

4.52(43)(63) 

4.45(43)(61) 


100.9 

133.3 

86.1 

116.8 

98.0 

128.8 

dof 

101 

113 

101 

113 

103 

115 

xVdof 

1.0 

1.2 

0.9 

1.0 

1.0 

1.1 


Table 5: Fit results for the dispersion relation matched to yPT at NLO. The renormalisation scale is /i = 770 MeV. 


The fit results for Lg are not in agreement with the determination of [40], 

10^ . Lg = 5.93(43). (152) 

Note that the influence of Lg on LJ, Lg and L 3 is minimal: if Lg is fixed to (152), we find 10 ^ ■L\ = 0.51(02)(06), 
10^ • = 0.89(05)(07) and 10^ • Lg = -2.82(10)(07). 

While the final results for the LECs do not differ significantly in the one-step and two-step strategies, a 
difference can be observed concerning the soft-pion theorems. If we use the two-step matching strategy, the 
soft-pion theorems are not automatically satisfied, but if they are imposed as a fitting constraint, they can 
be perfectly satisfied with only a slight increase of the In contrast, in the one-step strategy, where the 
subtraction constants have to fulfil the chiral constraints, the accuracy of the soft-pion theorems lies at ~ 4% 
and ~ 10% respectively. This does not change with the soft-pion constraints added to the fit, which only 
increases the x^ a bit- 

The influence of the isospin-breaking corrections of [28] is about the size of the statistical uncertainty in the 
case of Lg and Lg, while Lg is less sensitive to the isospin effects. 

A plot of the data points indicates that the two experiments NA48/2 and E865 are in agreement, which is 
confirmed by the fit results. We find it worthwhile to stress that this is only the case if the normalisation of the 
Fs data points of NA48/2 is determined including the s^-dependence (for the values in the 10 published bins, 
this requires the normalisation to be increased by 0.77%). If the published values are used, which are normalised 
neglecting the s^-dependence, a quite strong tension between the two experiments is observed, resulting in higher 
X^ values for combined fits. 

We note that the x^ in the dispersive treatment is clearly improved compared to the direct fit with xPT at 
NLO: in a fit to the one-dimensional data in appendix D, the x^ of a dispersive fit is 1.2 instead of 2.5 for the 
direct chiral fit (both with 27 dof). This is illustrated in figure 13: in contrast to a pure chiral treatment, the 
dispersion relation allows to describe the curvature of the S-wa,ve of the form factor F. We interpret this as the 
result of the resummation of final-state rescattering effects. Figure 14 shows the fitted P-waves of F and G. 
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5-wave of F 



Figure 13: Fit result for the S-wave of the form factor F. The dispersive description reproduces beautifully the 
curvature of the form factor. The (s,s^)-phase space is projected on the s-axis, the plotted lines correspond to splines 
through the (s, s^)-values of the two data sets. 


P-wave of F P-wave of G 




s/GeV^ s/GeV’ 

Figure 14: Fit results for the P-waves of the form factors F and G. The (s, s^)-phase space is again projected on the 
s-axis. 


5.2.2 Matching at 0{p^) 

We have seen that when using one-loop xPT, the dispersive treatment clearly exhibits its powers, and the 
advantage over a pure chiral treatment is evident: the dispersive representation is able to describe the energy- 
dependence of the form factors, hence the of the fit to the whole form factor data is much better. Due to 
the resummation of final-state rescattering effects, we expect the dispersive representation to capture the most 
important higher-order contributions and to render the determination of the LECs more robust. 

In combination with two-loop yPT, the treatment becomes more difficult. The matching equations at NNLO 
relate the subtraction constants to chiral expressions that contain the 0{p^) LECs C'*'. The largest obstacle in 
a chiral treatment at NNLO is the large number of poorly known Cl- In the dispersive treatment with NNLO 
matching, the same problem occurs. It turns out that the determination of the NLO LECs is still strongly 
affected by the choice of the Cl, a situation known from direct xPT fits [42, 37]. 

In order to alleviate this problem, we note that not all choices of the input Cl lead to a good convergence 
of the chiral expansion. In our dispersion relation, there appear nine subtraction constants, which are gauge- 
dependent quantities. Since the gauge transformation (74) is described by three parameters, we can find six 
gauge-invariant linear combinations of subtraction constants. For these linear combinations, we require a good 
chiral convergence. We obtain this by modifying the fitting procedure as follows. 

• We introduce 9 additional fitting parameters, corresponding to the contribution of the Cl to the subtrac¬ 
tion constants. 
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• We add to the nine observations of these parameters corresponding to the input values of the C’’ with 
a 50% tolerance for the linear combinations of the C^. 


• We add to the six observations of the total 0{p^) correction to the gauge-invariant linear combinations 
of subtraction constants. The observation is zero ±5.6% of the 0{p'^) contribution (5.6% corresponds to 
M^/iATTF^n 

With this setup, we are able to perform the NNLO matching with a reduced dependence on the input values of 
the Cl- In table 6, we present the matching results at NNLO, using the ‘preferred values’ of [37] as input for 
the Cf. 



NA48/2 

NA48/2 & E865 

NA48/2 

NA48/2 & E865 

10^ • L\ 
10 ^ • LI 
10 ^ • LI 

0.82(16)(09) 

0.71(10)(10) 

-3.10(40)(27) 

0.69(16)(08) 

0.63(09)(10) 

-2.63(39)(24) 

0.93(17)(04) 

1.11(17)(08) 

-3.96(49)(14) 

0.78(17)(03) 

0.97(17)(08) 

-3.38(48)(10) 

10 ^ • LI 

= 5.93 

= 5.93 

8.36(87) (48) 

8.05(86)(39) 


91.8 

123.9 

83.1 

115.3 

dof 

no 

122 

109 

121 

X^/dof 

0.8 

1.0 

0.8 

1.0 


Table 6: Fit results for the dispersion relation matched to xPT s-t NNLO. The renormalisation scale is p = 770 MeV. 
As in table 4, we use lattice input for LJ, L 5 , Lg and Lg [38, 39[, I /7 is the BE14 value [37[ and the input value for Lg is 
taken from [40[. 

The fit results with Lg taken as input are shown in the second and third column of table 6. Here, the 
corrections from NLO to NNLO matching for all three LECs are smaller than the corrections between NLO and 
NNLO observed in direct xPT fits. The larger uncertainties with respect to the NLO matching are explained 
by the additional fitting parameters for the C’’ contribution to the subtraction constants. If we take as input for 
the Cl the resonance estimate of [42[, we obtain {L^, L 2 , L^} = {0.65,0.26,-1.79} • 10“^. With the C} input 
taken from the chiral quark model [41[, we find {L},L 2 ;.^ 3 } = {0.49,0.65,-2.44} • 10“^. We prefer the BE14 
input values for the C}, because they lead to the best chiral convergence and the best x^ of the fit. 

The fit results change quite drastically if we include Lg in the fit. These fit results are shown in the fourth and 
fifth column of table 6. In the matching equations at NNLO, a stronger correlation between Lg and the other 
LECs is introduced due to their appearance in the s^-dependence. At present, alternative determinations of Lg 
are clearly more reliable than this one, and we therefore prefer here the fits with Lg taken as input. However, if 
the s^-dependence of the form factors will be measured in forthcoming experiments with even higher statistics, 
this could provide a new reliable way to determine Lg. 

5.3 Error Analysis 

In the following, we analyse the different sources of uncertainties in the determination of the LECs using the 
NLO and NNLO matching. Let us give once more the NLO and NNLO values for the LECs, obtained from the 
combined fits to the NA48/2 and E865 data, where Lg is taken as a fixed input [40[: 



NLO 

NNLO 

103.L}(/i) 

0.51(02)(06) 

0.69(16)(08) 

103.L^(/i) 

0.89(05) (07) 

0.63(09)(10) 

103.LS(/i) 

-2.82(10)(07) 

-2.63(39)(24) 


Table 7: Matching results for the LECs at NLO and NNLO. The scale is p, — 770 MeV. 


The first error indicates the statistical one, i.e. the error calculated in the linear fit of the parameters. This 
error is due to the uncertainty of the fitted data including isospin corrections. In the case of NNLO matching, 
it includes also the uncertainty introduced with the 50% tolerance of the C} contribution to the subtraction 
constants. The second error is due to the systematics of our approach. The corresponding statistical and 
systematic correlations are shown in tables 8 and 9. 
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total statistical 

form factor data 
isospin corr. 

C[ 

total systematic 

TTTT phases, low en. 
TTTT phases, high en. 
Ktt phases, low en. 
Ktt phases, high en. 

Ll 

LI 

Le 

Lj 

Ls 

L9 



NLO 


16.4 



NNLO 


Figure 15: Contributions to the uncertainty of Li in the 0{p^) and 0{p^) matching in units of 10 ® 


total statistical 

form factor data 
isospin corr. 

Cl 

total systematic 

TTTT phases, low en. 
TTTT phases, high en. 
Ktt phases, low en. 
Kn phases, high en. 



4.6 



NLO 


9.5 



Figure 16: Contributions to the uncertainty of L 2 in the 0{p'^) and 0{p^) matching in units of 10 ® 


total statistical 

form factor data 
isospin corr. 

Cl 

total systematic 

TTTT phases, low en. 
TTTT phases, high en. 
Ktv phases, low en. 
Kn phases, high en. 

LI 

Li 

LI 

Li 

Ls 

Li 



39.3 



NNLO 


Figure 17: Contributions to the uncertainty of L 3 in the 0{p‘^) and 0{p^) matching in units of 10 ® 
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stat. corr. 

r r 
^2 

TT 

-^3 

syst. corr. 

TT 

^2 

rr 

^3 

L\ 

0.49 

-0.72 

L\ 

-0.69 

-0.31 

LI 


-0.95 

LI 


-0.29 


Table 8 : Statistical and systematic correlations of the fitted LECs at NLO. 


stat. corr. 

TT 

^2 

rr 

^3 

syst. corr. 

rr 

^2 

T ^ 

^3 

LI 

0.31 

-0.32 

LI 

0.23 

-0.83 

LI 


-0.84 

L^ 


-0.70 


Table 9: Statistical and systematic correlations of the fitted LECs at NNLO. 

Figures 15, 16 and 17 show bar charts of the uncertainties of the LECs. The fractional uncertainties are 
summed in squares and determined as follows. 

• The uncertainty due to the Kn form factor data is the statistical uncertainty of a fit where no isospin 
corrections are included and the Cl contributions are fixed to the fitted values. 

• The uncertainty due to isospin corrections is the difference in squares of the statistical uncertainties of fits 
to data with and without isospin corrections, again with the Cl contributions fixed to the fitted values. 

• The uncertainty due to the Cf is the difference in squares of the statistical uncertainties of fits to isospin 
corrected data with the Cl contributions either fitted or fixed to the fitted values. 

• For the tttt phases [23, 24], we vary all the 28 parameters and sum the variations of the LECs in squares. 
In the bar charts, this is the uncertainty labelled by ‘tttt phases, low energy’. 

• The next fractional uncertainty is due to the high-energy behaviour of the tttt phases. We sum in squares 
the differences between the high-energy solutions explained in section 3.2.1. 

• The Kt: phases are simply varied between the centre and upper/lower limit of the error bands. This 
influence is labelled as ‘Ktt phases, low energy’. 

• The uncertainty due to the high-energy behaviour of the Ktt phases is estimated with the two solutions 
for each of the Ktt phases as explained in section 3.2.2. 

• The input LECs are varied by their uncertainties given in (151) and (152). 

• We have checked that the numerical uncertainties due to the discretisation, interpolation and numerical 
integration of the functions as well as the iteration procedure are completely negligible. 

We note that at NLO, the largest contribution to the systematic errors comes from the high-energy behaviour 
of the phase shifts, either from the tttt phases in the case of and L 2 or the Ktt phases in the case of L3. The 
uncertainties due to the low-energy parametrisation of the phases are small. The uncertainty due to the input 
LEC Lg is very small as well. 

At NNLO, the high-energy behaviour of the phases is again a large contribution to the uncertainty. Lg has 
now a large impact on the uncertainty of Lg. The additional LECs Lg, Llj and Lg have almost no influence 
on the uncertainty. The largest uncertainty is due to the fitted contribution of the Cl, which is part of the 
statistical uncertainty. 


6 Conclusion and Outlook 

We have presented a new dispersive treatment of Ki 4 decays, which provides a very accurate description of the 
hadronic form factors F and G. The dispersion relation is valid up to and including 0{p^) in the chiral counting. 
Furthermore, it provides a resummation of final-state tttt- and LfTr-rescattering effects, which we believe to be 
the most important contribution beyond 0 (p^). 

Our dispersion relation for K 14 is written in the form of an Omnes representation. It consists of a set of 
coupled integral equations. We have solved this system numerically with an iterative procedure. The solutions 
are parametrised by subtraction constants, which we have determined in a fit to data and by using the soft- 
pion theorem as well as chiral input. In contrast to a pure chiral description, the dispersion relation describes 
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perfectly the experimentally observed curvature of the S-wave of the form factor F, which we interpret as a 
result of significant TTTr-rescattering effects. This is yet another case in which high-precision data clearly call for 
effects which go even beyond NNLO in yPT. These effects only concern the momentum dependence of the form 
factors: we see no sign that quark mass dependence beyond NNLO is required by data. 

By using the matching equations to yPT we have extracted the values of the low-energy constants L\, L 2 
and L3. The correction from NLO to NNLO, when matching the chiral and dispersive representations and 
fitting the latter to the data are smaller than the corrections from NLO to NNLO observed in direct yPT fits. 
Constraints on the chiral convergence of the subtraction constants allow us to reduce the dependence on the 
input values for the C^. Still, the poorly known values of the are responsible for the larger uncertainties in 
the matching at NNLO. 

Our results for the LECs obtained by matching yPT at NLO are: 


L\ = 0.51(06) • 10-^ 

LI = 0.89(09) • 10-^ 

LI = -2.82(12) • 10-^ 

(153) 

whereas the matching at NNLO gives 

LI = 0.69(18) • 10-^ 

LI = 0.63(13) • 10-^ 

LI = -2.63(46) • 10"^ 

(154) 


The two-dimensional NA48/2 data set for the S'-wave of F, which shows both the s- as well as the s^-dependence, 
has allowed us to extract a value for Lg, which is roughly compatible with previous determinations. In accuracy, 
however, it cannot yet compete, as it reflects the low precision in the measurement of the s^-dependence of F. 
The determination of Lg is also quite strongly dependent on whether the matching is done at NLO or NNLO. 
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A Scalar Loop Functions 

We use the following conventions for the scalar one-loop functions in n dimensions: 

d^q 1 


Ao(m^) := V y 

:= V y 


( 27 r)” [<7^ — TO^] 
d^q 


1 


( 27 r)" [<7^ — ml] [(g + p)"^ — m|] 

These loop functions are UV-divergent. We define the renormalised loop functions in the MS scheme: 

Ao(m^) = —2w?\ + Ao(to^) -I- 0(4 — n), 

Bo{p'^,ml,ml) = -2A -I- Bo{p^F 0(4 - n), 

where 


(A.l) 


(A.2) 


(A.3) 


pL denotes the renormalisation scale. 
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The renormalised loop functions are given by [11] 

7^ 


o, m . ( m 

' — in 


Ao(m ) = - 


Boip'^ 


l) = - 


IGtt^ 

1 

IGtt^ 

1 






mt — mo 


2 

327r2 V ' \ ' A J \m? 


2 +-^ + - I In ' 


^ + vY - A^ 

“ p2 l^(p2_j,)2_ ^2 


where 


A 

E 


= mi - m2, 

2 I 2 

= rrii +1712, 


= A1/2( 


5 , 7711,1712 


)■ 


The renormalised two-point function fulfils a once-subtracted dispersion relation: 


Bo{s, m\,m\) = i3o(0, rn\, m\) -f - 


Imi3o(s^mf,m^) , , 


TT J^rm+m^r - S - ie)s' 


ds' 


where the imaginary part is given by 

Im_Bo(s, mi, m^) = 


2 _2^ ^ 1 A^/^(5,m^,mi) 

IGtt s 


and the value at s = 0 is 


Bo(0,mi,m2) = 


1 -m2ln(^) 


IhTT^ m\ — m; 

The first and second derivative at s = 0 are 

1 AE-2m2miln(^) 


Bp (0, mi, m2) = 


327r2 A3 

^ A(mf-I-10mfm|-I-m|) — hmim^Sln 


487r^ 


A5 


(A.4) 


(A.5) 


(A.6) 


(A.7) 


(A.8) 


(A.9) 


B Kinematics 

For each channel, the partial-wave expansion is performed in the corresponding rest frame, i.e. in the tttt centre- 
of-mass frame for the s-channel and in one of the Ktt centre-of-mass frames for the t- and rt-channel. Therefore, 
we work out explicitly the kinematics in the three different frames. 


B.l Legendre Polynomials and Spherical Harmonics 

For the partial-wave-expansions, we make use of several relations between spherical harmonics and Legendre 
polynomials. 

We use the addition theorem for the spherical harmonics and the relations between Legendre polynomials 
or derivatives of Legendre polynomials to spherical harmonics: 


Pi (cos 0') 

Pi'(cos 9”) 
P// (cos 6") sin 0" 


At: 

21 + 1 


^ Yr{e,Q)Yr*{e",<i7"), 


'1— — 1 


Att 


21 ' + 1 


Yl^(e",cl7") (for any 0"), 


t 7\ I ^ / 0-' + !)• vl* to" 

(-i)V2FTTV(F^^" 


(B.l) 

(B.2) 

(B.3) 
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Figure 18: Vectors and angles appearing in the addition theorem for spherical harmonics 


where Pl{z) := j^Pi{z)- The different angles are defined in figure 18. 

We can now easily derive the addition theorem for the Legendre polynomials: 

J dfPPi{cose')Pi,{cos9'') 

m——l 

= ^ J dn"Yr{e",r)Yi?ie",r) 




S,,,5„ 


4'7r / 47r 47r 

as well as the following relation: 

J dn"Pi{cos0')P[,{cose'') sin6i"e-*‘^" 

m=—l V ^ ' 

m——l V ^ ^ ^ 

"V" 

^rrtl 


Since the right-hand side is real, we conclude that 

f Ait 

/ dn"Pi(cos0')P//(cos0")sin0"cos(/)" = ^ii'w, -Ti^/(cos0) sinS, 

J 21 1 

J dfl"Pi (cos 0')Pl, (cos 9") sin 9" sin cp" = 0. 

B.2 Kinematics in the s-Channel 

In the TTTT centre-of-mass frame, the four-momenta of the different particles take the following values: 


(^y^M2. + fc2,fcj, 

qi = (^^/M^ + q^,q^ , 

Pl = { 

(-\jse + P,-k^ , 

to 

II 

to 

+ 

1 

P2 = ( 


^p) , 


(B.4) 


(B.5) 


(B.6) 


(B.7) 
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where <71 and 52 will be the momenta of intermediate pions. Note that we choose here the decay region is 
positive), but could have equally well chosen the scattering region. 

Inserting these expressions into s = {k — L)^ = [qi + = (pi +^2)^ gives the values of (f and fp'. We 

choose the directions of the three-vectors according to figure 19, i.e. the angles are defined as 9 := Z(—fc,pi), 
9' := Z(pi,gi), 9” := Z(-fc,gi). 



Figure 19: Vectors and angles in the s-channel centre-of-mass frame 
We end up with the following explicit expressions for the four-vectors: 


k = 


L = 


M'Y + s- Si \]^'j(s) 


2v^ 


M'Y - s- Si X^i (s) 


2^5 


2Vs 


2^5 


, 0,0 , 


, 0,0 , 


1/7 - M^cos9", J-- M2sin6l"cos(/)", J-- sin 6»" sin 


(B.8) 


92 = [ ^ - M2cos6»",-W)) - M2sin6»"cos())",-W)) - sin 6»" sin 


Pi = TV, 1/7 - M^cos9, M2sin6l,0 , 


P2 = 


s/I 

2 ’ 


- M2 cos9,-\ - - M2 sin0,0 , 


where Xki{s) := X{M‘Y,Si,s). Note that A]/J(s) has a square root branch cut in the s-plane between {Mk — 
■x/siY and {Mk + x/hY and changes sign when we continue it analytically to the scattering region. We will 
have to pay attention that we do not introduce this kinematic singularity into the partial-wave expansion. 

In order to express the s-channel scattering angle 9 with the Mandelstam variables, we calculate: 


t - u = {k - piY - {k - P 2 Y = + Pi - 2 kpi -k'^ - pi + 2 kp 2 

= 2k{p2-pi) = 2{k°{pl-pY -k - {P2-Pi)) 

= (-2y|^cos0) = -X][^^is)^l-^cos9 

= —2X(s)(T,r(s) COS0, 


(B.9) 


hence 


cos 9 = 


u — t 
2Xa^ 


(B.IO) 


where cr,r(s) := ^ 1 — 4M2 /s and Ar(s) 


hx]['^{s) as before. 
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B.3 Kinematics in the t-Channel 


In the t-channel, we are in the Ktt centre-of-mass frame and look at the t-channel scattering region: 


A:= WM|. + fc2,fc , qx = [ jMj. + 


-Pi= [yM^ + k‘^,-k] , q^=[jMl + qj^,-qK 


P2 = ( \lMl+f^,p2 
L = [ \/S£ -P2 I • 


(B.ll) 


Inserting these expressions into t = {k — pif = {qK + q-wY = (P2 + gives the values of P, ^ and We 
choose the directions of the three-vectors according to figure 20, i.e. the angles are defined as Ot := Z{—k,p 2 ), 
O't ■= AkAx), Q” '■= A-qK,P2)- 



Figure 20: Vectors and angles in the t-channel centre-of-mass frame 


We find the following results: 
k = 

Pi = 


qx = 


qir 


P2 = 


L = 


t +Ml- 

Ml 


or\G H, 
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5 
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9" cos A, - sin e” sin 4>” ] , 


t Mf. + ^XttA off ^xA) 

COS (7j , —^ sm (7j cos (p ^, 


2yft ’ 2v^ 

t-Si +Ml 

2^/1 ’ 2v^ ’ 

t +St- Ml 

2y/i ’ 2Vi 


2yft 


2Vi 


sm y" sin (/)" , 


0,0 , 


, 0,0 , 


(B.12) 


where \x-n{t) '■= A{M]p, Ml,t) and Xtnit) '■= A{st, Ml,t). Again, the square root of the first of these Kallen 
functions has in the t-plane a branch cut between {Mk — MtA and (Mk+Mt^)"^, the second between (Mn- — 
and (Mtt -I- Since we need the partial-wave expansion only in the scattering region t > {M^ + MA^, 

these branch cuts are not relevant. 

We calculate the t-channel scattering angle 0* as a function of the Mandelstam variables: 
s-u = {pi+ p 2 f‘ - {k- P 2 f =pI+pI + 2 pip 2 - A - pi + 2 kp 2 


— Ml — MI^ + 2p2 {k + pi) 


= Ml- Ml+ 2 


( t-st + Ml Ml-Ml 
y 2-\/t At 


2\/t \J t 



(B.13) 


hence 


cos Ot 


t(^S U) ^x-k^Itt 


(B.14) 
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B.4 Kinematics in the w-Channel 

The u-channel is completely analogous to the t-channel: 


k = 



-P2 = 




(B.15) 


Inserting these expressions into u = {k — p 2 )^ = {qK + qTr)"^ = {pi+L)^ gives the values of P, q^ and pf. We 
choose the directions of the three-vectors according to figure 21, i.e. the angles are defined as 6*„ := Z{—k,pi), 
'■= Ak,qK), OZ ■■= /L{-qK,Pi). 



Figure 21: Vectors and angles in the ti-channel centre-of-mass frame 
The results for the w-channel are then: 
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Let us calculate the u-channel scattering angle as a function of the Mandelstam variables: 

s-t={pi+ p 2 f - {k- piY =pI+pI+ 2 pip 2 -k"^ -p\ + 2kpi 
= Ml — Ml + 2pi(k -\- P 2 ) 

/ 

u 


= Ml-Ml+ 2 


-s, + Ml Ml-Ml Xli^u) X][liu) \ 
2^u ^ 2^ ^ cos„l. 
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(B.I 6 ) 
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C Omnes Solution to the Dispersion Relation 

C.l Solution for n = 3 Subtractions 


For n = 3 subtractions, the Omnes representation reads 
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Let us work out how to transform the Omnes representation (85) into the one with more subtractions (C.l). 
We start by subtracting all the dispersive integrals once more, using the relation 

1 _ 1 s 

s' — s s' (s' — s)s' (^-2) 

This generates nine additional subtraction constants: 
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To get rid of the subtraction constants in the i?-functions, we apply a gauge transformation (74). To this end, 
let us write the gauge transformation in the Omnes representation: 
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(C.4) 


Since the gauge transformation is a polynomial and has no discontinuity, the changes in the hat functions are 
given by SMq = —6Mo etc., which assures that the partial waves are unchanged. The shifts in the subtraction 
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constants are most easily found by comparing the Taylor expansion of (C.4) with (74): 
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where uj, to and Co are defined by applying subtractions to the Omnes functions: 
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In order to obtain the form (C.l), the subtraction constants in the i?-functions can now be removed with the 
gauge transformation 
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C.2 Hat Functions 


In the following, we provide the explicit expressions for the hat functions that appear in the Omnes solution to 
the dispersion relation. 
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sa-j^X 

~W 


4A4^ 


0-tX2 


- - (z^Ri)ty 

- yz^Ni)t, - {z^Ri)ty + ((2'‘Xi>t^ - {z*Ri)ty 

+ - {(««.).. -1'«.)..) si 


- ({z^Ni)t, - {z^Ri)ty 


2A/^7r — -5 + Sq 

4MJ 


+ ({2®Xi>t, - {z^Rih,) 


(7-n-X 

2Mj.' 
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lCT(f\-lAf\ ^Kn+t , ,, > {^iTT + t) ,,,> (Aktt +i) +t (So - 3t)) 

No{t) - {Mo)s, — - {zMo)s, - 


- {zMi)st 
{z^Mi)s 


(Aifyr i^iw(t) + ^l-K (A^tt + t)) + t i^lwii) + (^0 “ 3t) + t))) 






■ {Mi)s 


2t 


- (zMi), 






, t-Aj, 

' 6t 


,r N ./ON \AK^(i)(A/!,r+t) 
({zAf0>«t MzRo)ut) R,. 1/2.0 
DtA^^ (tj 


(A/f/r — t) (Ax/rA^^ + t (t - So)) (Ax,r A^^ + t (Sq - 3t)) 


+ ((-A^o)ttt “ 4(Ro)t, 

+ (^{Ni)ut — 4(-Ri}tit^ 

((zN ) A{zR ) ^ '^X 7 r(^) (A^tt + 0 (Axtt A^tt + t (/ — So)) (Axtt A^tt + ^ (^0 ~ 3t)) 

(Ax,r - t) A]/^(t)A^^/^(/) (Ax/r A/.^ + *2)' 


24/3 M|. 


12/3M|. 


//2 at-\ / 1 / 2 p\ ( Axtt ^) Axtt (t)\e^{t) >^Kir{t) (A^^ + /) (Ax,rA^.r + t^) 

+ ((. /Vl}„. - 4(. j -j 

(A/.+t) 


24/3 M|. 


_ ^ — 4(_Ri} ^ ^Axtt (3A^^ + Sq + /) + A|^^ (A^^ — 2/) + 3/2 (Sq — /) 


+ (^{zNl)ut - A{zRl)nt^ 


- ((z‘^Ni)ut -A(z^Ri)u,) 


2At‘^Mj^ 

x]ll{t) (Ae^ + /) (Ak^ {Ae^ - 2/) + / (Sq - /)) 
2At^Mj,\li\t) 

(Ax/r+3/) A]y/^(/)A^^/^(/) 

24/2 m2. 

Ax/r)*) (A^.n- + /) 

24/2 m2. 
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iVi(t) 


(Mo), 


3Ml,{Ae^+t) , , 3M^(AK.+t) oa^4 


9Mf^ (A£„ +1) 


StXUt) 8tA,^(t) 


3M2.(A^^+t)(A^,A^^ + t(So-3t)) 

• (Mi)s4 

8t (^) 

3Mj, {Ak„ + t) (Ajy^A^^ + t (So - 3f)) _ (A^^ + t) 

U^x]llit)xli\t) 8t2Ai/"(t) 

9M|.(A^,+t)(A^,A^^ + t(So-3t)) 3M2,(AK^+t)‘ 




8t^Xe^{t) 


4t2 


9^1:+ 


■ (2 “l)st --T79-^ (™l)st --77^-1- (^™l) 

4tXe„{t) 

f^i.-ry i^llT + i) f/AT \ -1/p' '\ i^i-T + t) 

^ 4tA,,(t) MRo)u,) 


3Mj^ {Aktt - 3t) 

2tx]^:(t)xli\t) 


UXi^it) 


+ {{zNo) ut - 4:(zRo)u 
+ {{Ni)u, -4{Ki>„,) 


M^jt-AKn) 

■^tx]{l{t)xli\t) 


+ ({z^No)ut -4{z^Ro)ut) 


3Mj^ (A(^ + 1) 
4tXe^it) 


{^£n + i) {^Ktt^t/tt -\- t (t — So)) {Ak-rAit^ + t (Sq — 3t)) 


+ ((2iVi)„^ -4(2 Ri)„^) 


16t3A^^(t) 

(Ajf-jr — t) +t{t — So)) (AjCtt + t (Sp — 3t)) 

(*) 

^KttV) i^t-TT + i) (AifTr^-Z-n- + 




8t^xli\t) 

3 (A^.^ + 1) (Aj^rAir -\-1 (t — So)) {AkrAir + t (So — 3t)) 

16t^XeR{t) 

^Kirit) {^e-iT + t) {Akr — t) {AkrAir + t^) 


16t3 


4t3 


I //.3i\r' 1/.3p \ X [ i^K-n- R)^KRi^)^lL (*) . 3A^^(t) (A^.^ + t) (Ax-n-A^^ + 

+ ((• '''■>«. - “I' «.)..)(-^3--^-- 

- {(•■‘Wi)., - 40*R,+ *> 


16*3 


-4 {(S,).. - KS.),.) «1i^-+<)(‘‘ + ”^r.-':o<-^r..‘A-) 

I (/~\T\ -IZ-P ' \ f(^iTT + t) 

+ {{zN,U -MzR,U) 


+ ({22iVi)„, -4(22fii>„^) 


(tAxn (3A^7r + So + i) + (A^tt ~ 2t) + 3t^ (Sq — i)) 

8tn](lit)xli\t) 

3iW^ (A^tt + i) (Axtt (A^tt ~ 2t) + t (So — t)) 


16t2A£7r(^) 


(Axtt + 3t) 
8^2 


^.3.-.' 1X3B1 N3A4^A;/Ai)(A*R+i) 

- (0 iv.)., -40 R.)..)--, 


(C.12) 


64 



~ o, 3iW^ 9 3 “t“ t t^o — 3t)) 


U 


- ((1 - z^)zM,)., 3A^/^(t)A,^/^(t) ^ 


4t2 

3M| 

2t 


- (((1 - z^)NoU - 4((1 - z^)Ro)u,) 

+ (^((l-z^)Ni)„, -4((1-22)Ki}„^) 


i^K-rr^tw +t{t— Eq)) (Axtt^^tt + * (^0 — 3i)) 




- (((1 - z^)zN,U - 4((1 - z^)zR^U,) W(Ax.A,.+t^) 

+ (((1 - z^)z^mu, - 4({1 - z^)z^R,U,) 

_i_ /"/('I ‘^\AT\ A /("I 2\ p \ \ ^ ‘^^^Ktv ~ ^0^ ~ 

+ (((1 - 2 )^l}ut - 4((1 - 2 )Rl}ut) - 


+ (((1 - z‘^)zNi}ut - 4((1 - z‘^)zRi}^ 


8^2 

8t2 


(C.13) 


Ro(«) = (Mo>.„ ^ + - (.A4o>.„ + («(So - 3u) + Ak.A,^} 


4u 


Auxyj^iu) 


Su^M% 


+ {zMi)s 
- (z^Mi) 


(“)4^^(“)(“ +, {u +Ae„) (uC£o - 3u) + Aii^Ae^) \]l^{u) 


8«2m2, 


Su^M^xyjiu) 


2,^\ + Alf7r(«) , > ^U-Aktt , , , (« + («) 




™ 4uA^; (tj) 

I _ (2(.iVo>*„ + 

6« V ^ GuXyyiu) 


(^Ktt - t*) (“ (“ — So) + AjfTrAe^) (u (So — 3u) + AKn^in) 


+ ( 2 (Afo}t„ + {RoK 
- (2(Afi}t„ +(i?i}t„) 

I (n/^AT \ I /»p \ A / (^Ktt — u) xy^{u)xy {u) [v? + AkttAit,) 

+ {2{zNr)t^ + (zR,),^) - 


24«3m|. 


^ (ti + AtT^) (u (u — So) + Ap^ A^^) (u (So — 3n) + Ap^ A^^) aJ^^(«) 
24u3M^xj^^(u) 

i/o/2Ar\ I / 2 D \ \ f ~ ^Ktt) ('*^)-^^7r ('^) i'^ 4“ ^^tt) ("^ + ^KTri"^) A 

+ {2{z 7Vi>*„ + (. j -- j 

+ (2o»iv.)..+(.»«.),.) 

+ (2(7Vi}t + (_Ri}t ) 3 (So -m) + Ak ^ (u + So + 3Ae^) u + A2,^ (A^^ - 2u) 


24u‘^Mj^ 


fo/ AT \ I / D \ \ ( ^ SqW + 2Ap7r'*^ ^KtT^Itt} 

+ e20«.).,. + OR,)..) (-24.«««A;f(„,- 


+ 2(^27Vi}t„ +(22Ri}t^ 


(3m + Ay^) A]/^(tt)A^)^^(tt) 
24u'^M]^ 

("U H- A/^7j-('U) 

24«2m|. ’ 
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Rliu) = {Mo)s^ + (zMo) ——K 


- {Mi)s, 

+ (zMi) 

+ {z^Mi)s^ 


8u\e^{u) 

3M|, {u + Ai^) (u (So - 3u) + Ak^Ai^) 


(u -\- Aktt) /^ 2 ! 17 -. \ ^ 

- 1/2, x\l/2/ ., ~ -^o)s 


8uXi^{u) 


IG-u^ A^tt (w) 

3M]^\]^^{u) (A{^ + It) ^ 3M|, (Ajf^ + u) {Ak-k^i-k + m (Sq - 3n)) ^ 

16«2Aj/^(it) 8u2A^^(n)Aj/^(ti) 

9M2, (A(.^+«) (AkitA^.^+«(S o - 3u)) 3M2,(AK7 r+«) 


IG-u^ A^tt (-u) 


8ii^ 


9M|. (« + A,,) Aj^/2(n) 


- (.*M.).. - (A},).. _ ,.A},) “"K 




3Mj. (Axir - 3n) 


16M2A^/^(n) 


8nA^^(«) 




f.2tZr\ S-^xC^ + A^ir) /r,/,^r \ I /P' -^K (“ + A<,r) 

(n/^AT \ I /-p \ A (Axir — t*) (o1^2at\ i /~2p\ N 3M^ (« + A^^) 

(n + A^^n-) (it (u — So) + Aktt^Itv) (u (So — 3tl) + Axir A^^n-) 


16u^\eT,{u) 


- (2(iVi)t„ +(iJi)t„)' 

+ (2(^7Vi)t„ + (2/?i}t„) [ W (A^. + u){Ak^A,^ + 

_ (Aktt — (Aktt + ^ (^ - So)) {Aktv^Itv + ^ (So - 3n)) \ 
3 (AfT^ + It) {AK-TT^tir +u{u- So)) (AxttA^^ + tj (So - 3ti)) 


+ {2{z^Ni)t^ + {z^Ri)t^'^ 


16u^\e^{u) 

>'KTr{u) (Ai^ + u) (Axir - m) (AxttA^^ + «2) 


16u^ 


4u^ 


+ (2(.3iVi}*„ + (^3iJi},„) 

_ 3A^^(u) (A^^ + u) (Axtt A<;,r + ^2) 

, A^x (“ + A^ir) (u (So - It) + Axir (A^,r “ 2u)) 
16u^Xe„(u) 


(2(iVi>t„ + (Ri>t„ 

(" 2 ( 27 ^ > + (zK > ') f A 42 , (uAxvr (3A<;^ + So + u) + A2^^ (Afa - 2u) + 3tt^ (Sp - u)) 

^ " 1 ‘aiJ 8n2Aj/^(«)A]/^«) 


+ 2(z^Ni)t^ + (z^Ri)t, 


^K-^KnM (A^ir + «) 

16u^xj^^(u) 

3M2, (A^^ + m) (-Axu A^^ + 2uAxir + - Spu) 

16^2 A^tt (-u) 




(Axir +3n)\ 

8«2 J 

(A^vr+^) 

16n2A]/^(u) 
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where 


and 


^ f \ //I 2^ 71\ 2^ 71^ \ ^ (^0 ^w)) 

Ri(u) = {{1- z )Mo)su —;-((1-2 - 


4u 


8ii^ 




+ ((1 - 3 ^ 

+ (2((1 - ^2)7Vo),^ + ((1 _ z^)Ro)t^) ^ 

- (2((1 - z^)Ni)t^ + ((1 - z^)Ri)t^'j {AKn^i^+ui u - So)) + «(So - 3n)) 




/ 0 O \ ("w) (^iCTr-^^TT + '^^) 

+ ( 2((1 - z^)zN,)t^ + ((1 - z^)zR^)t^) \uiMl -^ 

- (2((1 - z^)z^N^)t^ + {(1 - z^)z^R^)t^) ^KAn)\lA^) 

+ (2((1 - z^)Ni)t^ + ((1 - z2)Ri>t„) 

- ( 2((1 - z^)zfR)t^ + ((1 - z^)zRi)t^'^ 


Aktt (A^tt - 2m) + u (So — u) 


8n2 

8«2 


{z-X),^ :=1 J"^z-X{t{s,z))dz, 

{z'^X)s, ■=\J z'^x{s{t,z))dz, 

(z"X)„, := 1 j\^X{uit,z))dz, 

{z"‘X)s^ ■=\ J z"-X{s{u,z))dz, 

(z-X),^ :=1 J\-Xit{u,z))dz, 

t{s, z) = 1 (So - S - 2Xa„z), 

s{t,z) = ^ ^So - t + 1 (^zX][l{t)Xyj^{t) - AK^A(>yj'^ , 
u{t,z) = 1 ^So-t - 1 (^z xyi{t)xyy{t) - AK^Ai^j^, 

siu,z) = ^ ^So - ■“+ 1 , 

tiu,z) = 1 ^So-u- 1 (^z Xyi{u)X]y{u) - Ak^Ai^'^^ . 


We recall the abbreviations 

Aktt = Xlf^ — Af^, Ai^ = Si — M^, So = + 2Af^ + si, 

PL =-{M]( — s — si), X = -X^^‘^{M](, s, si), cttt = -^1 — ^ ^ . 
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D Isospin-Breaking Corrected Data Input 

In this appendix, we list the isospin-corrected data sets on the form factors that we use for the fits of the 
dispersion relation. These are the NA48/2 [6, 9] and E865 data sets [7, 8], corrected for isospin-breaking mass 
effects and (in the case of NA48/2) the additional radiative effects that were calculated in [28]. 

More detailed explanations can be found in section 4.1. 

D.l One-Dimensional NA48/2 and E865 Data Sets 


v^/MeV 

V^/MeV 

Fg 

Fp 

Gp 

Gp 

286.06 

92.61 

5.6941(85)(185) 

-0.181(67)(15) 

5.035(257)(66) 

4.317(74)(20) 

295.95 

92.01 

5.7878(90)(170) 

-0.324(62)(34) 

5.168(142)(84) 

4.404(53)(32) 

304.88 

91.51 

5.8410(89)(171) 

-0.209(60)(33) 

4.924(108)(59) 

4.532(46)(26) 

313.48 

90.65 

5.8905(91)(171) 

-0.156(58)(32) 

4.879( 91)(51) 

4.627(41)(24) 

322.02 

88.32 

5.9275(90)(166) 

-0.366(55)(41) 

5.227( 80)(58) 

4.692(38)(29) 

330.80 

85.59 

5.9557(93)(168) 

-0.383(54)(40) 

5.265( 73)(56) 

4.748(35)(28) 

340.17 

81.02 

5.9915(92)(166) 

-0.218(55)(46) 

5.036( 68)(59) 

4.762(34)(31) 

350.94 

76.16 

6.0161(92)(163) 

-0.302(54)(35) 

5.246( 62)(37) 

4.889(34)(21) 

364.57 

69.80 

6.0351(91)(162) 

-0.309(54)(33) 

5.338( 57)(31) 

5.000(35)(20) 

389.95 

58.96 

6.1155(93)(224) 

-0.264(59)(35) 

5.400( 55)(34) 

5.144(36)(22) 


Table 10: NA48/2 data [6, 9], corrected by additional radiative and isospin-breaking mass effects [28]. The uncertainties 
of the isospin corrections (without the higher order estimate) are added in quadrature to the systematic error. The fully 
correlated error of the normalisation increases from 0.62% to 0.70%. The normalisation of 77, is increased by 0.77% to 
take the dependence on si into account. 


V^/MeV 

yi^/MeV 

Fg 

Gp 

287.6 

106.8 

5.781(13)(42) 

4.702(89)(40) 

299.5 

105.7 

5.825(14)(48) 

4.693(62)(37) 

311.2 

103.8 

5.914(14)(56) 

4.771(54)(41) 

324.0 

101.1 

5.974(16)(62) 

4.999(51)(56) 

339.6 

96.3 

6.097(17)(63) 

5.002(49)(57) 

370.0 

84.6 

6.151(20)(41) 

5.104(50)(42) 


Table 11: E865 data [7, 8 ], corrected by isospin-breaking mass effects [28]. The uncertainties of the isospin corrections 
(without the higher order estimate) are added in quadrature to the systematic error. The fully correlated error of the 
normalisation is 1.2%. 


D.2 Two-Dimensional NA48/2 Data Set 

For the fits of the dispersion relation to data, we do not use the above NA48/2 data set on Fg consisting of 10 
bins, but the two-dimensional data set, which was recently published as an addendum to [9[. Here, we list the 
isospin-corrected values of Fg that we use as input in our fits. The values and uncertainties, shown in table 12, 
are constructed as follows. 

• With the number of data and Monte Carlo events for each 2D bin [9], we compute the relative values and 
statistical uncertainty of the relative values of Fg. 

• We fix the normalisation by requiring fg = 5.705 in a parametric fit of the form (112). 

• Unfortunately, systematic errors are not available for the 2D data set. We guess a systematic error by 
assuming that the ratio of systematic and statistical error does not depend on si. 

• We apply isospin corrections due to photonic and mass-difference effects [28]. The uncertainty from the 
mass effects is added in quadrature to the systematic error. 
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Fs 

1 

2 

3 

4 

5 

6 

1 

5.641(51)(55) 

5.628(30)(35) 

5.700(24)(30) 

5.687(22)(28) 

5.744(21)(27) 

5.707(22) (28) 

2 

5.753(51)(30) 

5.716(30)(22) 

5.747(24) (20) 

5.777(22)(20) 

5.807(22)(20) 

5.833(23) (20) 

3 

5.785(52) (33) 

5.775(31)(24) 

5.831(25)(21) 

5.837(22)(21) 

5.841(22)(20) 

5.875(23)(21) 

4 

5.908(52)(33) 

5.809(31)(23) 

5.877(25)(21) 

5.874(22)(20) 

5.910(22)(20) 

5.894(23) (20) 

5 

5.910(52) (24) 

5.903(31)(19) 

5.891(25)(18) 

5.909(22)(18) 

5.924(21)(18) 

5.961(22)(18) 

6 

5.831(51)(29) 

5.925(30)(22) 

5.912(24)(20) 

5.919(22)(19) 

5.971(21)(19) 

6.032(22)(19) 

7 

6.031(50)(25) 

5.927(29)(20) 

5.941(23)(18) 

5.970(21)(18) 

6.024(20)(18) 

6.045(22)(18) 

8 

6.026(47)(21) 

5.976(28)(18) 

5.990(22)(17) 

6.020(20)(17) 

6.024(20)(17) 

6.067(23)(17) 

9 

6.023(44) (22) 

5.987(26)(18) 

6.037(20)(17) 

6.059(19)(17) 

6.044(20)(17) 

6.077(25)(18) 

10 

6.163(38)(67) 

6.128(22)(41) 

6.107(18)(34) 

6.120(18)(34) 

6.139(23)(42) 

6.130(45)(79) 


7 

8 

9 

10 



1 

5.703(25) (30) 

5.721( 30)(35) 

5.717(43)(47) 

5.709( 82)(86) 



2 

5.817(26)(21) 

5.828( 31)(23) 

5.872(43)(27) 

5.929( 84)(45) 



3 

5.843(26) (22) 

5.934( 31)(24) 

5.911(43)(29) 

5.923(117)(68) 



4 

5.905(26)(21) 

5.957( 31)(24) 

6.111(50)(32) 




5 

6.004(25)(18) 

5.942( 33)(20) 

6.074(70)(29) 




6 

6.025(26)(20) 

6.009( 38)(24) 





7 

6.042(28) (20) 

6.124( 54)(26) 





8 

6.086(33)(19) 

6.024(122)(39) 





9 

6.058(54) (25) 







Table 12: Values of Fs for the two-dimensional data set of NA48/2 [9] including isospin-breaking corrections [28]. The 
fully correlated error of the normalisation of 0.70% has to be treated separately. 


E Matching Equations 

E.l Subtraction Constants at 0{p‘^) in xPT 

In the following expressions for the subtraction constants at NLO, we have used the Gell-Mann-Okubo (GMO) 
formula = (4M|- — M^)/3 to simplify the analytic expressions considerably. This introduces an error only 
at NNLO. In practise, we use the physical rj mass and not the GMO relation. We do not show the analytic 
expressions for this case because they are much larger. 
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'*^O.NLO 


Mk 


d,2L\Mf. + SLIM% 




■ 273M|.M2 - 258M|.M4 - 


-49M^) 


- In 


+ ln 


- In 


3847r2A|,^(4M|. - M2) 

Mj-iSMl- - 24M|.M2 + 21Mj^M* - 7M|) 

M2,(38M|. - 6M|.M2 - 9M2.m4 + 3M6) 


/i2 


M| 


m2 


M^ 


647r2A|^^ 

(4M2, -M2)2(5M2, +M2) 

128^2^3^^ 


0 

-i,Nno = ^ 


- 8L^M|, 

M2,(79M|. - 2M2,m2 + 7M4) 


-In 


■ In 


- In 


/i- 

/i2 


384^2a^^ 


/A7|\ 5M2,m4(3M2, - m2 ) 




128^2a3^^ 


M2,\ m6.(43M2. - 21M2) 


647r2A3 

i\ TT 


M2,(4M2, - M2)2(5M2, + m2) 


128^2^3^^ 


(E.2) 


(E.3) 


- 0 

"^l.NLO = 


V2F„ 


1 + rj^ ^ — SL 2 {M'^ + 2M2 + se) — 4 L^{M]^ + m2) + 4L5M2 + 2 LqSi 


16M6 - ZM^Ml ■ 


■3M2,m4 


- In 


■ In 


- In 


/,2 


M^ 


■2M| ^ (M2,+M2)2 

647r2A2^^ 


967r2A2^^ 

m2 (3M2 - m2 ) (m4 - 8M2 m2 - 5M4 ) M4 (3M2 - m2 ) 

-S<! 


128^2^3^^ 


M2,(37M6. - 35M^M2 - 17M2.M4 + 3M6) 


64^2a3^^ 


M2 \ 68 M|. + 7M^M2 - 2M2,m4 - M^ \ 


1287r2A2^ 

i\ TT 


327r2A3„^ 

i\ TT 


M^(M2,-3M2) 

327r2A3,^ 


(E.4) 


"’■l.NLO — 


Mk 

V2F3 


SL^M^ 

M2,(M^ + M4) 


32n^Al^ 


- In 


- In 


m2 \ MUM% - + 12M2 M4 - 4M6) 


967r2A3^^ 


M2,^ M2,(7M|. - 21M|M2 + 3M2,m4 - m 6 )\ 


1927r2A3 

A. TT 


(E.5) 


1 _ Mk 

"■O.NLO “ 


v^F3 


- 24L5M|- - 6 L 3 MI- 


M2,(16613M®. - 2179M^M2 + 29Mj^M^ + mM^) 


- In 


■ In 


- In 


M2 

//2 


20487r2A3^^ 


M2\ 3M2 M4(37M4 - 80M2 m2 + 20M4) 


5127r2At- 

A TT 

m2, ^ 3M|,(-4840M4, + 1216M2.M2 + 83M4) 


3M2,(4M2, - m2)(304M|. - 6M|.M2 - M^) 


128^2^3^^ 


(E.6) 
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E.2 Matching at NNLO 

E.2.1 Decomposition of the Two-Loop Result 

E.2.1.1 NLO Contribution 

We have already decomposed the NLO contributions. We apply a gauge transformation to convert the expres¬ 
sions to the second gauge and evaluate the result numerically: 


m. 


0,NLO 

0,L 


m 


m 


m 


m 


n, 


n, 


M, 


K 


1, NLO _ 

2, NLO 
^0,L 

0,NLO 

1,L 

l,NLO 

l.L 

0,NLO 
^1,L 


V2F^ 

_ Mk 
V2F^ 
= 0 , 

^ Mk 
V2F^ 
= 0 , 

Mk 


l.NLO 

1, L 

2, NLO 

1,L 

1. NLO 
0,L 

2, NLO 
0,L 


V2F^ 

Mk 


V2F^ 

= 0 , 

Mk 


V2F^ 

0,N 

^1,L 


( - 0.1466 • -f 0.4953 • lO^L^ -f 0.0872 • lO^L^ -f 0.0733 • lO^L^ 


0.0092 • lO^L 


3 Tr 


0.0573 • lO^LS^ 
9 Ml 


(o.9173 • lO^L^ -f 0.2293 • lO^L^), 


( 


- 0.2293- lO^L^,, 


( - 0.2660 • lO^L^ - 0.1238 • lO^L^ -f 0.0092 • lO^Lg 
- (0.2293 • lO^L^ - 0.0573 • lO^Lg)^), 

( 0 . 


2293 • lO^L^ ,, 


( - 0.6880 • lO^L^ - 0.1720 • lO^L^), 


O.NLO ~l,NLO n 
= n,'r =n{j^ = 0, 


(E.7) 
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m 

m 

m 

m 

m 


0,NLO 

o,R 


l,NLO 

0,R 


2,NLO 

0,R 


0,NLO 

1,R 


l,NLO 

1,R 


m 

rh 

rh 


0,NLO 

1,R 


l,NLO 

1,R 


2,NLO 

1,R 


n 


l,NLO 

0,R 


n 

n 

n 


2,NLO 

0,R 


0,NLO 

1,R 


l,NLO 

IM 


Mk 

V2F^ 

Mk 

Mk 

Mk 

V2F^ 

Mk 

V2F^ 

Mk 

V2F^ 

Mk 

V2F^ 

Mk 

V2F^ 

Mk 

V2F^ 

Mk 

V2F^ 

Mk 

V2F^ 

Mk 

V2F^ 


(0.1393 + 0.0444^+0.0256-|-), 


0.3413- 0.0512 


si 


(0.4080), 


-0.0916- 0.0512 


(0.0512), 


si 

Ml 


-0.0902 - 0.0595^- 0.0256^,^, 


0.1545 + 0.0512 


S£ 

Ml)' 


(0.0137), 


-0.1376- 0.0384 


St 

Ml)' 


(-0.0796), 
(0.0384), 
(-0.0282). 


(E.8) 
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E.2.1.2 NNLO LECs 

First, we consider the contribution of the NNLO LECs, the . We decompose this contribution into the form 
of the polynomial part in (69): 

^ ^ (4^^ (Cl - 2^3^ - 2CI + 2Q + 4Q 

" " + 2C[o + 8Cri - 4C[2 - 8C[3 + 202^2 + 402^3 " 2C'J4) 

- Ml Ml (4C'[ + 64C'2^ + 56C3^ + 34Q - SCg^ + 40C^ 

+ 640^ + 2ACI - 16C'[o - 48C'[i - 8C'[2 + 112C'C3 
+ 16C'[4 - 80C[5 - 64C'[7 + 80^2 - + 160^5 

+ 16C2"6 + 32C2"g + 64C^o + 32(7^6 

I /^r I o/^r /^r , /^r \ 

"T ^66 ^^67 ~ ^69 ~ ^88 ^90/ 

+ Ml{- 24C'[ - 12802^ - 3203^ - 18C; - 320^ - 2406^ - 64Cf 
+ 8Q + 8C'[o + 16C'[i - 80C'r2 - 80C'C3 + 32C'r4 + 8C[5 
+ 128C'[g - 48C'[7 + 160^2 + 3202^3 + 3202^6 + 128C'^8 

r//or I o/^T /^r i 

~ ^1^66 "T ^^67 ~ ^69 ~ ^88 ^90 

+ (mI (4C[ - %Cl - QCl - 8C'r2 - 32C'r3 

O^r _ 1 _ _|_ I O/^T /^T Q/^'^ \ 

~ 0^63 ~ "T ^66 "T ^^67 ^^^69 ~ ^88 ~ '^^90/ 

- Ml (12C[ + 64C'2^ + 72CI + lOC^ - 480^ - 

- 32CI, + 8Cl, + 4^4 + 4^5 + 9Qg + 2Cl, 

+ 16^8 + 3^9 + 8^3 + 16^4 + 3^8 + Qo)) 

- 2sl (SCI + 2^1 + Qe + 2^7 - Qg + Qs - Cg^o)) , 

^ ]r, (sMl (Cl + ACl + 4^3^ + 2CI + 4Q + 4^7^ + 2Q 

" " +2C[2+8Cr3-2C2^3 + C2^5) 

+ 2MlMl (2ACI + <d&Cl + 32CI - &Cl + SCI 

+ SCI + 16C^ - 16C'[3 - 16(7^3 + 802^5 

I /^r I /^r /^r i /^r \ 

~c ^66 ~c ^^67 ^69 ^88 “*"^90; 

+ 2s(Ml (4C[ + 16CI + 16C'3^ - 2Cl 

+ 3Qe + 2C^,, + 4Cg^8 + Cg-g + C;, - C;,)^, 

^ 16M4 (-d- 4d - CJ + Q), 

^0;NNLO ^ (mi (4^3^ + 2CI - 4C[o - 16C[4 + 12C[2 + 32C[3 

— ^^22 ~ ^^23 ~c ^^63 ^66 ~c ^67 ~ ^69 ~ ^^83 ~c ^90) 

+ Ml (I6C3" + I6Q - 4C[o - 8C[i + 4C[2 + 16C[3 

- 10^2^2 - 802^3 + 6C2^5 - 2^3 + ^Cl, + 2C8^3 - Qg) 

+ s,( 12C'3^ + 2CI + Qg + Cl, - Cl, + Cls - Qo)). 

= ^^2Mi( - 16^3^ - 6^ + Cg^g - Cg^7 - Cg^g - + C^9o). 
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m 


0,NNLO 

1,C 


^ SCI + lOQ - 8Q - 8Q - 8C[o - 16CI, + 16CI, 

^ ^ ^ + lec^a + 8C[5 + 16C'[7 - 16(7^2 - 32C'^3 

I /^r I o/^r /^r i /^r 

"T ^66 ^^67 ~ ^69 ~ ^88 ^90/ 

+ M^mI- ( - 12C'[ + 2ACI + 34C'I - SC^ - 24^6^ - 8CJ 

- 8C'[o - 32C'ri + 16C'r2 + 166^3 + 16C'[4 + l&Cl^ 

- 36(7^2 - 32(7^3 + 20(7^5 + leC^g - 32^2^9 

+ 4^3 + Qe - 6^7 - Qg - 4^3 + Qs + Cg^o) 

- 2M|. (2C[ - 2(7^ + ACl + 8(76^ + 12C'[2 + 16C[3 + 2(7^2 

+ 4C2^3 + 4CJ4 + 2^3 + Qe + Qt - Qg - + Cqo) 

+ (mI (- 4C'[ + 40(73^ + 2&CI - 8C'[o - lec'j'i + 8C'[2 + lec 

- 12^2^2 - 16^2^3 + 4C2^5 - 4^3 ” 4^4 ” 4^5 + ^66 

— QCqi — SCgg — 4Cg3 + Cgg — TCgg) 

- Ml (4C'[ + 8C3^ - QCl + 8C[o + 32C[i - 16C[2 - 32C 

+ 4C2"2 + 8C2"3 + 4^3 + 8^4 + Qg + 2^7 + 3Q, 


13 

r 

69 


+ 4Cg"3 + Cg"g + G 


88 


-'90 


m 


l,NNLO 

1,C 


1 

V2F^F^ 


2,NNLO _ Mk 1 


n.c 


l,NNLO 


'-o,c 


y2F^ 

1 

V2F^F^ 


+ s^,{8C;-2C^,l 
2Ml (-Ml {4C; + 8Cl - 4C[o - leCi^i + 4C'[2 + 32C'[3 - 6C;. 

- AC;, + AC;, - 2 C'g ^3 - 2 Cg ^7 + 2 Cg ^3 - Cg^g) 

- Ml (lecj + lOQ - 4C[o - 8Cl 

+ 12C[2 + 16C[3 - 6C2"2 - 8C2"3 + 2C2"5 

I I /^r _i_ ^ 

-h ^003 -h Oge — ZOg 7 — L/gg — ZOg 3 -h ^goj 

- s,(12C3^ + 2Cl + Cg^g - 2C;, - Cl, + c;, - Cg^o)), 

o /i/f4 f Q/^'^ _i_ _i_ _i_ \ 

ZiVJ^ 1^503 -h ZO 4 -h L/gg — Og 7 — Ogg — Ogg -\- Oggj, 

8Ml ( - 2Ml {8Cl - ACl + 2Cl + 4Q 

+ 2C[o + 8Cri - 2C[2 - 8C[3 + 2C2"2 + 4C2"3) 

- (^Ml (16C[ + 8Cl - 18CI + 8Cl + 80^ 

+ 8C[o + 16C[i - 16C[2 - leCi^g + 16(7^2 + 32(7^3 

0/^4’ I I /-ir \ 

~ ^66 ~ ^^67 ^69 "T ^88 ~ ^90/ 

-\s, {ACl - SC^, - 6Cl + Qg + 2^7 + 3 Q 9 


-^8 + ^ 0 )), 


2,NNLO _ _ Mk 1 


'-o,c 


0,NNLO 

^1,C 


.l,NNLO 

^l.C 


V2F^Fl 

Mk 1 -3M^ 


12 M^(Cr + CJ-CJ), 


iC 


V2F^Fl 
Mk 1 
V2F^Fl 


( 1I /^r 0/^4" I /^r 1 /^r 

l^iOOg -h DO 4 — L/gg — ZL/g 7 -h L/gg -h L/gg — L/ggj, 


8Ml {8Cl + 2Cl + Qg + 2Cl, - Cl, - Cl, + Cl,). 


(E.IO) 


Unfortunately, a lot of NNLO LECs enter the polynomial. In total, there appear 24 linearly independent 
combinations of the Cf. 

If we use the resonance estimate of [42], we obtain the following values for the NNLO counterterm contri- 


74 



bution: 


O.NNLO 

Mk 

^0,reso 

V2F^ 

l,NNLO 

Mk 

”^0,reso ~ 

V2F^ 

2.NNLO 

Mk 

^0,reso 

V2F^ 

O.NNLO 

Mk 

^l.reso 

V2F^ 

l.NNLO 

Mk 

"^l.reso ~ 

V2F^ 

~ O.NNLO 

Mk 

”^l.reso ~ 

V2F^ 

~ l.NNLO 

Mk 

^l.reso 

V2F^ 

~ 2.NNLO 

Mk 

"^l.reso ~ 

V2F^ 

l.NNLO _ 

Mk 

^O.reso 

V2F^ 

2.NNLO _ 

Mk 

^O.reso 

V2F^ 

O.NNLO _ 

Mk 

^l.reso 

V2F^ 

-l.NNLO _ 

Mk 

^l.reso 

V2F^ 


-0.1546 - 0.1716^ + 0.0316^ 
Ml Ml 


0.1747- 0.0316- 


(E.ll) 


(0.0310), 

1 0.1657- 0.0316^ 

V M 

(-0.0104), 
^-0.0900-0.0135- 
^-0.1712-0.0316- 
(0.1805), 

1 0.1502 - 0.0237^ 

V M 

(-0.0233), 
(-0.0078), 

(0.2707). 


Alternatively, if we use the ‘preferred values’ of the BE14 fit [37] (complemented with Cgg — Cgg = —55 • 10 ® 
[40] and the remaining LECs that appear in the s^-dependence set to zero), we obtain the following values for 
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the NNLO counterterm contribution: 


0,NNLO 

"^0,BE14 


l,NNLO 

"^0,BE14 


2,NNLO 
^ 0 ,BE 14 


0,NNLO 

^ 1 ,BE 14 


l,NNLO 

"^ 1 ,BE 14 


~ 0,NNLO 
^ 1 ,BE 14 


~ l,NNLO 
"^ 13 E 14 


~ 2,NNLO 
"^ 13 E 14 


l,NNLO 

^0,BE14 


2,NNLO 
^0,BE14 


0,NNLO 

’^ 13 E 14 


~l,NNLO 

^ 1 ,BE 14 


Mk 

V 2 F^ 

Mk 

V 2 F^ 

Mk 

V 2 F^ 

Mk 

V2F^ 

Mk 

V 2 F^ 

Mk 

V2F^ 

Mk 

V2F^ 

Mk 

V2F^ 

Mk 

V2F^ 

Mk 

V2F^ 

Mk 

V2F^ 

Mk 

V2F^ 


f-0.4108 - 0.1823^ - 0.0033^^ , 

V Ml Ml)' 


0.7959 + 0.0986 


(-0.1709) 


0.2627 + 0.0296 


St 

Ml 


St 

Ml 


(-0.1709), 

f0.0356 + 0.1050^ + 0.0263^^ , 

V Ml Ml)' 


-0.2942 - 0.0296 


(0.1841), 


St 

Ml)' 


0.3505 + 0.0296 


St 

Ml)' 


(0.0099), 
(-0.1282), 
(0.2761). 


(E.12) 


E.2.1.3 Vertex Integrals 

Let us study in more detail the contribution of the vertex integrals. They can be decomposed into functions of 
one Mandelstam variable according to 


ciNNLO 

tv 


{s,t,u) = Flsfi°{s,st) 


+ 


u — t _nnlo 


NNLO 


/ \ , i-ilNlNEI 

(s,s^) -\-Fvt,o 


(i, s^) + 


Ml 


niNNLO/. \ , 77NNLO/ \ 
■FyT,l {t^S^)-\-Fyu [U^Sl) 


^NNLO 

\jy 


(s, t^u) = G 


NNLO 

VS 


(s, Si) + G 


NNLO 

VT,0 




+ 


s — u 


^NNLO 

^VT,l 


{t, St) + G 


NNLO 

VU 


{u,st). 


(E.13) 


The M-channel vertex integrals fulfil In the following, we treat them numerically. The 

contribution to i?o is obtained by subtracting the constant, linear and quadratic terms: 


(u, s,) = St) - St), 


pNNLO 

^VU 


{u,st)=Fll 


NNLO 


{Q,st) + uFll^^'{Q,st) + -yF, 


2 pNNLO" 


VU 


( 0 , 5 ^), 


(E.14) 


where ' stands for the derivative with respect to the first argument (u). The polynomial Py^'"® has to be 
lumped into the overall polynomial and finally reshuffled into the subtraction constants. Numerically, we find 


-0.2521 - 0.0130^) i + 0.0569 

As we have checked again numerically, the polynomial-subtracted u-channel contribution of the vertex integrals 
fulfils the dispersion relation 


P 


NNLO 


VU 


{u,st) 


Mk 

V 2 F^ 


0.4008 + 0.0119 


St 

Ml 


Ml 


TTT ■ (E- 15 ) 


Rq (■“> = 



ImRVu'.s,) 

y — u — 


(E.16) 
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Next, we consider the s-channel vertex integrals: apart from a polynomial, they belong to either Mq, Mi or 
Ml. Again, we subtract the first few terms of the Taylor expansion: 


mY (s, se) = Si) - (s, Si), 

MY(s, Si) = Si) - PG,ys°(s, Si), 

Ppy^sois, Si) = -Fys,o°(0: Si) + sJ^vs,o°'(0, Si) + ^s^A^s,o°"(0> Si), 

Ppy'fiis, Si) = aNnl°(0, Si) + sF^^^.^'iO, Si), 

j’Sf (»,«) = Gsy°(o,«) + + is"G!;s“"(o, «>. 


(E,17) 


We find numerically 


pNNLO 

^F.VSO 


{s,Si) 


pNNLO 

^F,VS 1 


{s,Si) 


pNNLO 

^G,VS 


{s,Si) 


Mk 

V2F^ 

Mk 

V2F^ 

Mk 

V2F^ 


0.2663 + 0.0992^ + ( -1.7763 - 0.0450 

M]c 


Sj 


0.0029 + 0.0006^ + 0.0006 

Mk 


Ml 


-0.3197 - 0.0727^ + [ 0.1457 + 0.0163 
Ml 


Si 

'Ml 


Ml 


- 0.5385 


Ml 


- 0.0003 


Ml 


Ml 


(E.18) 


A numerical check shows that the polynomial-subtracted s-channel contributions of the vertex integrals fulfil 
the dispersion relations 


Mlis,si) 


Ml (s. Si) 


Ml{s,si) 


si r°° lmMl{s',Si) 
TT Jso {s' - s- ie)s'^ 
s2 /•“ ImMi^(s',s^) 
^ Jso {s' - S- ie)s'‘^ 
si /■“ ImMi^(s',s,) 
^ Jso {s' - s - i€)s'^ 


(E.19) 


Finally, we consider the t-channel, which is a bit more intricate: the reason is that not all linear and quadratic 
terms of a simple Taylor expansion in t belong to the subtraction polynomial. The f-channel contributions can 
be written as 


pNNLO 

FvT.O 


{t,Si) 


nNNLO 

^VT,1 


{t,Si) 


^NNLO 

^VT,0 


{t,Sl) 


^NNLO 

^VT,l 


{t,Si) 


, 2AKn^en,.Vu “2 Ak-^ - St ~y 1 dV"/. , dNNLO 


-Nr{t,si)-- 


3 2M 


g-No (t,s,) + -- 

2 t 

3AT_ 

-Ik {t, Si) - {t, Si) + ^ KhKK (b Si) - {t, Si) + Si), 


2 K {t, Si) + 7 ’F,VT°(t) Si), 


2 -Ni {t,Si) + -Ro {t,Si) + Pp^VTo{t,Si), 
K G 


2Ml 


K{t,Si) + PlK[{t,Si), 


SMI 


(E.20) 


where y¥o second order and Pa^vri order polynomials. The Taylor expansion 

of Nq starts with a cubic term, the one of with a quadratic and the one of NY with a linear term. 
Numerically, we find 


pNNLO 

^F,VT 1 


{t,si) = -piiKi 


{t, Si) 


Mk 

Kf^ 


( 0.0044-k 0.0002^ 

V Ml 


0.0003 



(E.21) 


and also identify the linear and the quadratic term of the Taylor expansion of NY . In the sum 


pNNLO 

^VT,0 


{t. Si) -f G 


NNLO 

VTfi 


{t, Si) 


At 


T^K (^. + PFyTo{t^ Si) + PlyK*’ 


NNLO , 


(E.22) 
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we can easily separate from the sum of the polynomials. After having identified particular the 

quadratic term of its Taylor expansion), we can also separate the difference of the polynomials using 


s,) - s,) = it, Si) + i 


'iNNLO 


4 


Mi 

jNNLO 


Nr{t,se)-- 


2 Aktt — t 


3 Mi 


NYit,Si) 


+ {t, Si) + Si) — Pa^Toit, Si) 


Numerically, we find 

oNNLO/. . \ . Mk 


-){t,Si) 


oNNLO /, ^ >. 
^G,VT0[t: Si) 


V2F^ 

Mk 

■n/STtt 


-0.6831 - 0.1136 


-0.0055 - 0.0146 


St 

Mi 

St 

Mi 


-0.0013 


- 0.0006 


Mi 

St 

Mi 


+ 0.2841 - 0.0006 


St 


Ml Mi 


^ +0.0190 ^ 


Again, the following dispersion relations can be checked numerically: 

+3 /•- lmN^(t',Si) 


.rv, s F f In 

^oit,Si) = -l - 


+ 0.0131 + 0.0095 


dt', 


St 

Ml 


Ml 


t t 

- 0.0356^ 
Ml Ml 


/3 


(t' — t — ie)t‘ 
ImNlit', Si) , 




Nl{t,si) = - 

TT 


io 


IrnAj (t',Si) 

(t' — t — ie)tP 


dt'. 


Reshuffling the polynomial contributions into the subtraction constants leads to 


m, 


O.NNLO 

OM 


m, 


l.NNLO 


OM 


Mk 

V2F^ 

Mk 

V2F^ 


0.0705 + 0.0667yy - 0.0539^^ 


Mi 


Mi 


-1.7841 + 0.0648 


,2.NNLO ^ ^ (-0.5954) 
V 2+7r 


'■OM 


m 


O.NNLO 

1 .+ 


m 


n, 


n, 


l.NNLO 

^i.y 


. O.NNLO 
^i.y 

. l.NNLO 
^i.y 

2.NNLO 

1 .+ 

1. NNLO 
0,V 

2. NNLO 
0,V 

O.NNLO 

^i.y 

-l.NNLO 

^i.y 


Mk 

V2F^ 

Mk 

V2F^ 

Mk 

V2F^ 

Mk 

V2F^ 

Mk 

V2F^ 

Mk 

V2F^ 

Mk 

V2F^ 

Mk 

V2F^ 

Mk 

V2F^ 


Si 

Ml 


Si 

Ml 


-0.2658 + 0.0969 


(-0.1132) 


-0.1310 - 0.2580^ + 0.0435^ 

Mi Mi 


0.2570 - 0.0849 
(0.0572), 

-0.2587 + 0.0519 
(0.0898), 

(-0.0849), 

(0.0729). 


Si 

Ml 


Si 

Ml 


(E.23) 


(E.24) 


(E.25) 


(E.26) 
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E.2.1.4 Remaining Two-Loop Integrals 

Next, we consider the remaining two-loop parts, It is easy to decompose them into functions of one 

Mandelstam variable: 


77 NNLO 


{s,t,u) 


(^NNLO 


{s,t,u) 




NNLO/ ^ 
PS 




NNLO/ 


S — U 


PT,0 ,^^2 ^PT,1 (tjSe) 


Ml 


'rNNLO/„ „ \ I piNNLO 


FP {s,t,u), 


s,) + s,) + s,) 

+ Gl^^^iu,Si) + PE!l^^is,t,u), 


(E.27) 


where and are second order polynomials. The remaining steps are analogous to the case of the 

vertex integrals. Again, we apply subtractions to the different functions: 


Mlis, se) = s,) - Pll'fis, s^), 


M[ (s. Si) = Si) - Pg.ps° (S) si), 

RI {u, Si) = Si) - P^pp° (m. Si) 


_ piNNLO 
— ^PU 


iu,Si)-PB%^ 


{u,Si), 


(E.28) 


where 


pNNLO 

^F,PS 


(s, Si) 


pNNLO 

^G,PS 


(s, Si) 


pNNLO 
^F,PU 


(u, Si) 


pNNLO 

^G,PU 


(u, Si) 


G?g“(0.«) + »G!;g“'(tl, „) + i»»G!;g“"(0,«), 

G?y°(o,„)+„G!!y<>'(o, „) + i«=G!!y<>"(o, »,>. 


(E.29) 


Numerically, we find 


oNNLO 

F,PS 


{s,se) 


oNNLO/ \ 

Pg,ps (s,se) 


Mk 

V2F^ 

Mk 

V2F^ 


-0.1660 + 0.0002^ + 0.0007^ -t ( 1.1629 -t 0.0343 

Ml Ml 


nNNLO/ \ nNNLO/ \ 

Pf,pu (U,Sl) — Pg,pu (u,si) 


0.0609-t 0.0118 


Mk 

V2F^ 


St 

Ml 


+ -0.0514- 0.0058 


St 


St 

Ml) Ml 

2 


^+0.7815^ 


Ml Ml 


® -0.0007 ^ 


Ml 


St 


0.0585 - 0.0089^ -t 0.0165 -t 0.0069 

Ml 


St 


Ml Ml 


“ - 0.0442 “ 


Ml 


The t-channel contributions can be written as 


t-iNNLO / / \ ^ AjP t-i- \ I ^ ^^KlV^^tlV TijP ii \ 2 ^^Ktt 3t ^jP/i \ I I T)P 14- \ I t^NNLO / j 8 
FpT,o (b St) — g -^0 (^) St) + - —— Fi (t, St) — - — 2 M 2 —Pf,pto(^i st), 

J^NNLO/, 8 2t ,,rP/, 8 1 77)NNL0 7, 8 

PpT.l (tjSt)— g^2 -^1 P^St) + Pp^PTlitjSt), 

^NNLO/, 8 atP t-i 8 2 AKn^tn ,.jP,, •, , 2 Ap'.n- -|-t ,r,p , , I p, , r,NNLO /, 8 

GpTfi (t, St) — — -Nq (t,st) — - ——Ni (t, St) -I- g 2M ^—~ 3"^° PG,PTo\i, st), 


Gpp] 7 °(i, St) — “3^2 (^’ St) -f Pc^p^Tiit, St), 


^K 
^NNLO, 


(E.30) 


(E.31) 


where PppT®, Cs pto second order and Ppp^i -^G p¥i order polynomials. Numerically, we find 


F,PT1> G,PT1 


Pl^P^^iit, Si) = -P^^p^Hit, Si) « (^-0.0010 - 0.0001 


Ml 


and also identify the linear and the quadratic term of the Taylor expansion of . In the sum 


(E.32) 


4t 


J^PT.o (i, Si) = Si) + Si) + PSTtW^ 


(E.33) 
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we can separate Nf from the polynomials. We obtain the difference of the polynomials with 

4 AjPu \ ^ ~ t -\jP 14 . \ 


’ O O IVlj^ O ivij^ 


+ ^^0 ~ ^G,p¥o(^5 


(E.34) 


and find numerically 


P 


.NNLO 


F,PT0 


(t, Si) « 1 0.2047 + 0.0339^ 

V2P V 


_0.1781 +0.0019j|-)j^-0.0211j| 


P, 


f-0.0662 + 0.0085-®^ 


G.PTOO-) ~ ^ -u.uuuz, -r 


p: 


Si! 


0,0757-0,0054^J ^+0,0251^ 


Finally, the additional polynomials are given by 




iNNLO 

F.P 

iNNLO 


Reshuffling all polynomial contributions into the subtraction constants leads to 


(E.35) 


(E.36) 


m 


O.NNLO 

0,P 

1. NNLO 

\,P 

2, NNLO 

\,P 

O.NNLO 

'^l,P 

l.NNLO 

■^l.P 


m 


O.NNLO 

l.P 


m 


l.NNLO 

l.P 


m 


2,NNLO 
l.P 


n 


n 


1. NNLO 
0,P 

2, NNLO 
'^O.P 

O.NNLO 

l.P 


n 


l.NNLO 

l.P 


Mk 

V 2 F^ 

Mk 

V 2 F^ 

Mk 

V 2 F^ 

Mk 

V 2 F^ 

Mk 

V 2 F^ 

Mk 

V2F^ 

Mk 

V2F^ 

Mk 

V2F^ 

Mk 

V2F^ 

Mk 

V2F^ 

Mk 

V2F^ 

Mk 

V2F^ 


0.3349 - 0.0426^ + 0.0429^ 


Ml 


1.1922 - 0.0523 


(0.8257) 


Sj 

Ml 


-0.0083 - 0.0797 


(0.0884), 


Si 

Ml 


0.0771 + 0.0016^ - 0.0367^ 

Ml Ml 


-0.0030 + 0.0757 
(-0.0449), 
-0.2217-0.0478 
(-0.0693), 

(0.0662), 

(-0.0633). 


Si 

Ml 


Si 

Ml 


(E.37) 


E.2.1.5 NNLO One-Loop Integrals 

The last NNLO piece that we have to decompose is the part containing the U). Similar to the two-loop parts, it 
can be easily decomposed into functions of one variables. Since this contribution contains only one-loop integrals, 
we can express it in terms of Aq and Bq functions, which can be treated analytically. After decomposing the 
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NNLO one-loop part according to 


t, u) = 

+ Fi!^^0(u,se) + FP^^^^(s,t,u), 
t, u) = G2r°(^, + G2 ^.o°(^, 5^) + 


(E.38) 




NNLO/ 


the polynomial contribution is found in analogy to the two-loop part. Reshuffling the polynomial gives very long 
expressions for the subtraction constants. We perform a Taylor expansion in S£ and evaluate the expressions 
numerically, using the physical masses and p. = 770 MeV: 

O.NNLO _ Mr 


‘0,L 


\/2Fn 


(0.0243 -I- 0.0155 ■ 10®Lg) • 10®LJ -|- (0.3528 - 0.0523 ■ 10®Lg) ■ 10®LJ 

+ (0.0831 - 0.0092 ■ 10®LJ) ■ 10®LJ -I- (0.0400 -I- 0.0350 ■ 10®LJ - 0.0020 ■ 10®LJ) ■ 10®LJ 
+ (0.0066 -I- 0.0048 ■ 10®Lg) ■ 10®LJ - (0.0012 -I- 0.0699 ■ 10®LJ -I- 0.0087 ■ 10®LJ) ■ 10®Lg 
+ 0.0213 ■ 10®Lg -I- (0.0100 - 0.0027 ■ 10®Lg - 0.0003 ■ 10®Lg) ■ 10®Lg 
+ ^ (0.0213 ■ 10®Lg - 0.0161 ■ 10®Lg -I 0.0230 ■ 10®Lg -I- 0.0139 ■ 10®Lg -I 0.0018 ■ 10®Lg 
- 0.0017 ■ 10®Lg - 0.0008 ■ 10®Lg -I (0.0229 - 0.0060 ■ 10®Lg) ■ 10®Lg) 


+ ( - 0.0053 ■ 10®Lg - 0.0029 ■ 10®Lg - 0.0029 ■ 10®Lg -I- 0.0065 ■ 10®Lg 

-I 0.0010 ■ 10®Lg - 0.0012 ■ 10®Lg - 0.0006 ■ 10®Lg -I- 0.0025 ■ 10®Lg) j , 

^i.NNLO ^ f _ (0.1644 + 0.0968 ■ 10®Lg) • 10®Lg - 0.2921 ■ 10®Lg - (0.1665 -I- 0.0242 ■ 10®Lg) ■ 10®Lg 

- 0.0353 ■ 10®Lg -I- 0.0049 ■ 10®Lg + 0.0185 ■ 10®Lg - 0.0033 ■ 10®Lg + 0.0076 ■ 10®Lg 

+ (o.0138 ■ 10®Lg - 0.0575 ■ 10®Lg - 0.0087 ■ 10®Lg - 0.0130 ■ 10®Lg 

- 0.0020 ■ 10®Lg -I- 0.0024 ■ 10®Lg -I 0.0012 ■ 10®Lg -I- 0.0196 ■ 10®Lg) j , 

.^2,NNLO ^ (^0.3345 • 10®Lg -I 0.2734 ■ 10®Lg -I- 0.1618 ■ 10®Lg -I- 0.0863 ■ 10®Lg 

+ 0.0096 ■ 10®Lg -I- 0.0067 ■ 10®Lg - 0.0003 ■ 10®Lg -I- 0.0032 ■ 10®Lg j , 

.^O.NNLO ^ (" _ g J203 ■ 10®Lg -I- (-0.2247 -I- 0.0242 ■ 10®Lg) • 10®Lg - 0.0727 ■ 10®Lg 

- 0.0241 ■ 10®Lg - 0.0046 ■ 10®Lg + 0.0078 ■ 10®Lg -I- 0.0039 ■ 10®Lg 

+ (0.0121 ■ 10®Lg -I- 0.0044 ■ 10®Lg -I 0.0063 ■ 10®Lg - 0.0130 ■ 10®Lg 

- 0.0020 ■ 10®Lg -I- 0.0024 ■ 10®Lg -I 0.0012 ■ 10®Lg - 0.0053 ■ 10®Lg) j , 


l.NNLO _ Fix 

“ V2F^ 


- 0.0198 ■ 10®Lg - 0.0059 ■ 10®Lg - 0.0056 ■ 10®Lg -I- 0.0130 ■ 10®Lg 


+ 0.0020 ■ 10®Lg - 0.0024 ■ 10®Lg - 0.0012 ■ 10®Lg 


(E.39) 
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- 0,NNLO 

T = ~i =— 

V2F, 


0.0440 ■ 10^1,5; + (-0.1488 + 0.0281 ■ lO^Lg) ■ lO^LJ + (-0.0140 + 0.0131 ■ 10®L5) ■ 10®Lg 

+ (0.0001 + 0.0044 ■ 10®Lg) ■ 10®L4 + (-0.0033 + 0.0048 ■ 10®Lg) ■ 10®Lg 

+ (0.0186 - 0.0087 ■ 10®Lg) ■ 10®Lg - 0.0135 ■ 10®L7 + (0.0026 - 0.0003 ■ 10®Lg) ■ 10®Lg 

+ [ - 0.0957 ■ 10®LJ - (0.2423 - 0.0242 ■ 10®Lg) ■ 10®LJ - 0.0520 • 10®Lg - 0.0134 ■ 10®LJ 

Mr ^ 

- 0.0033 ■ 10®LJ + 0.0067 ■ 10®Lg + 0.0033 ■ 10®Lg + (0.0098 - 0.0060 ■ 10®Lg) ■ 10®Lg) 

g2 

+ (^0.0057 ■ 10®LJ + 0.0010 ■ 10®LJ + 0.0029 ■ 10®Lg - 0.0065 ■ 10®LJ 

- 0.0010 • 10®Lg + 0.0012 ■ 10®Lg + 0.0006 ■ 10®Lg - 0.0034 ■ 10®Lg) j , 

^i.NNLO ^ _Mj£_ ( 0.0987 ■ 10® 1,5; + (0.2328 - 0.0242 • 10®Lg) • 10®L^ + 0.0581 ■ 10®LJ + 0.0213 ■ 10®L); 
v2i^7r \ 

+ 0.0062 ■ 10®Lg - 0.0078 ■ 10®Lg - 0.0039 ■ 10®Lg 

+ ( - 0.0138 ■ 10®Lg - 0.0029 ■ 10®Lg - 0.0026 ■ 10®Lg + 0.0130 ■ 10®Lg 

+ 0.0020 ■ 10®Lg - 0.0024 ■ 10®Lg - 0.0012 ■ 10®Lg + 0.0089 ■ 10®Lg) j , 

^2.NNLO ^ / _ 0.0070 ■ 10®Lg + 0.0114 ■ 10®Lg - 0.0097 ■ 10®Lg - 0.0044 ■ 10®Lg 

v2i^7I- \ 

+ 0.0001 ■ 10®Lg + 0.0012 ■ 10®Lg + 0.0006 ■ 10®Lg^, 

^i.NNLO ^ Jp£_ ^ _ Q Q^gg . (_0.4712 + 0.0726 ■ 10®Lg) • 10®Lg + (-0.1097 + 0.0181 ■ 10®Lg) ■ 10®Lg 

V 2i^7r \ 

- 0.0262 ■ 10®Lg - 0.0049 ■ 10®Lg + 0.0075 ■ 10®Lg - 0.0117 ■ 10®Lg - 0.0021 ■ 10®Lg 
+ ^ (0.0095 ■ 10®Lg + 0.0035 ■ 10®Lg + 0.0040 ■ 10®Lg - 0.0098 ■ 10®Lg 

- 0.0015 ■ 10®Lg + 0.0018 ■ 10®Lg + 0.0009 ■ 10®Lg - 0.0079 ■ 10®Lg) j , 


2,NNLO _ 


0,NNLO _ Mk 


\/2LV 


- 0.0003 ■ 10®Lg - 0.0010 ■ 10®Lg + 0.0101 ■ 10®Lg - 0.0007 ■ 10®Lg 


- 0.0022 ■ 10®Lg + 0.0003 ■ 10®Lg - 0.0010 ■ 10®Lg - 0.0004 ■ 10®Lg 


- 0.0125 • 10®Lg - 0.0051 ■ 10®Lg - 0.0069 ■ 10®Lg + 0.0098 ■ 10®Lg 


+ 0.0015 ■ 10®Lg - 0.0018 ■ 10®Lg - 0.0009 ■ 10®Lg ), 


-l.NNLO 
^1,L 


\/2LV 


0.0059 ■ 10®Lg + 0.0096 ■ 10®Lg + 0.0051 ■ 10®Lg - 0.0072 ■ 10®Lg 
V 

- 0.0012 ■ 10®Lg + 0.0018 ■ 10®Lg + 0.0009 ■ 10®Lg j. 

Note that there are no quadratic terms in Lg, Lg or Lg. 


(E.40) 


E.2.2 Chiral Expansion of the Omnes Representation 

In order to derive the NNLO chiral expansion of the Omnes representation (C.l), we first expand the Omnes 
function chirally: 


n 


NNLO 


f X ^ S _ S® 


<5nlo(s') 


{s' — s — ie)s‘ 


/3 


ds' 


+ 


1 


Mi 


s 

TT 


<5lo(s') 


{s' — s — ie)s'^ 


,ds' 


(E.41) 


where the subtraction terms to and Q are defined in (148). 

In the quadratic term of the expansion, only the LO phase enters and therefore only two subtractions are 
needed. The NLO expansion of the modulus of the inverse Omnes function is given by 


1 _ s 

|^NL0(5)| = — 


5^ho{s') 

{s' — s — ie)s''^ 


,ds'. 


(E.42) 
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Therefore, the NNLO chiral expansion of the argument of the dispersive integrals reads 


M(s) sin(5(s) 


|f}(s)| 


= M^°(s)5nlo(s) + M”^^(s)(5lo(s) 


NLO/ 


NNLO 




This leads to 


,^NNLO^„',_ Mo/^, ,0 ^ , ,-0 I ^0,NLo(^') 


ds' + - ( uJq -I - 


'5o,lo(^0 j ^ 
as 


+ I Aa"“ + 6““ — + c"“ —^ I w°— , ^ A” '^o,Lo(«') 
+ I ^“NLO + "nLO + %LO j,^4 j I <^0 ^2_ + 




_L ^1^0 —_U r^O _ 

^ “NNLO ^ '^NNLO j^2 ^ ^NNLO pj-4 

jijNNLO f _ f^M, ,Mi ^ ^ f ,. .1 ^ I 


(s'-s-is)s'^ 


‘^l,Lo(^0 
(s' — s — ie)s'‘ 


tcIs' 


,2 ,3 r“Mr°(s')^},LoL'),, 


-I- o“i -I- ^_L r“i A_1_ :_ / , , / 1.1.W, / ,„/ 

T t^NNT.O ''nNT,0 1, rS ^ f4 i, rA ^ I , . ■ \ ,3 ’ 


-NNLO V -NNLO V -nnLO ^^4 ^ ^ J ^ _ .^y. 


,VNNLO/„',_ Mi/^, ,1 « , ,-,l f” ^l,NLo(^0 .„/ 1 / ,1 

■ “"O P A4 + T y,„ (s'-s-*.y" 2 P 


^ r '^LLoL') , , 


, All - 111- , - I ‘^l,Lo('®0 

NLO ^ “NLO + %LO j^4 ) I “1 M| + ^ y,„ (s' - S - *e)s' 


Ml , 7 Ml ^ . ^Mi 


+ ^^i^T.n + 


NTPT + - 


^ds' 


I ^Mi . 7 Ml 
^ “nnt.o ^ “ 


-NNLO -r -NNLO + %NLO ^^4 T -NNLO Mf, TT J(s' - S - ie)s'‘^ 


■ + 


s3 ,4 r“Mr°(s')-5;,Lo(*'), , 


^ H- 


ds , 

^oihoi^') 


YNNLO/jN ( jNo ^ I gJVo p ^0°(*')<5oyo(*') A / 1/2 t f , -u.LOV/ ^ , 

^^0 W-^^Vo]^^+%LO]^^ + 7/^ (t'-t-ie)t'^ '"PP ^ + ^y*„ (LLTTpp'^* 

+ 6^0 p ^0^°(oClo(O 

■r ^NNT.n J 1 rO "t" ^NNT.n ji ,«-4 "r / /., . . x . ,fl / - • . , .o tli- 


-NNLO - -NNLO ^ ^ _ -y 


r N^°(t')6y,loin ,« <' t- 


(t'-t- u)n 

*„ (t'-t-nn v'-'’ Ml' (t"-t'-ie)t"^ ^) ’ 

oo \ t r^r°(t'pLo(0,. 


i+p 


AtNNLOu\ - ..Ni f,,l/2 ^ ‘^l,Lo(^') , „iVi , - / L X 2 

^ W-^nlo(^^i M|"^7ri (t'-t-ie)t'^ {t'-1 - ie)t^ ’ 

YNNCO/j-, _ y^iVi * , p ^l'°(*')<5iyo(*') A /_l/2 t I P ‘^LLo(*') 

(i) - [b..ojy + -1 j Mi + L i it'-t-ie)n^^ 

, i , r ^i°it')^i.Loin /“ A"^°(t')AfLo(<') 

+ f'NNLo^il + ^A^ (t'-t-nn + (t'-t-oi'" 


<2 /■” ^J'°(t')po(0 L 1/2 t' t'" 


-, 1/2_ 

jt„ (t'-t-nn V Afi 2t„ 

/*! r^o°(oCo(opp3/2 t r ^o.Lo(o 


jt, (t" - 1 '- it)t"^ j 


-dt a;. 


/ta (t'-t- ie)n ) Ml TT Ao (t'-t- ie)t' 

*3 iir(t')<A'Lo(*')r ^r°(i')^of/o(t') 

+ ' (t--t-ie)n (t'-t-nn 


dt' 


Jto {i- - i- tt)f 'i Jto 

^ /°° Ar (<')Co(t') ^ _,3/2 t' , t'^ 


R^^^°(t) = 0, 

^NNLO(j) ^ 


(f — t — it)f 


/3 


Ml 


Co(i") 


(t" - t' - ie)f" 


idt" 1 dt' 


where we use the following notation for the contributions to the subtraction constants: 


(E.43) 


(E.44) 


OnLO = OLO + AonlOi 
ONNLO = OLO + AonLO + AonNLO- 

Remember that and apQ are non-zero after the gauge transformation. 


(E.45) 
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We further define 


uMa — 

"nlo “■ 


,0 


M 




^NNLO 


=: -Wn 


Mk 

V 2 F^ 


+ b 


NNLO> 


f^NLO 


lMi 

"NNLO 


-^;^ + bS,- 


'V 2 F, 


NLO’ 


Mk 


_l_ 


which allows the simplifications 


_ Mk ( f 0 1 02^ .s^r < 


'5o,NLOv“ / , / , 

as + 


(^') 

- ie)s'^ 


1 IS 

2 


'5o,lo(*0 

(s' - s - ie)s'^ 


^ds' 


+ I Aa"" + 6"“ — + c"“ —^ — + — n '^o.Lo(^') ^ A 


Mr''°(s) = 


Mo , rMo ^ 
"■NKLO ^ ^NNLO j^2 


Ml ,Mi ^ 

^NLO ^ ^NLO j^2 


Ml I rMi 

“NNLO ^NNLO j^2 


■NNL0^4 +'*NNLOjlf6, + ^ 

^ ^ <^;,L0(''') 

^ JT y,^ (s' _ s - je)s'^ 

3 rx, A>NrO|'„'Ui, 

J sn 


f 

J Sr\ 


,ds’ 


[s' - s - it)s'‘^ 


Mf^^°(s) = 


Mk 

F 2 F„ 


-NNL0^4 T ^ 

„2 „3 r°° (5J 

(s' 


■Mr°(s>;,Lo («')., 

-ds , 


/ OO si 

ds' + 

„ (s'-‘'-’A"'" 


^ 1 ^“NLO ^ "NLO j^2 


"NLO ^4 


(s' — s — ie)s'^ 

^^iis'+ V- 
ie)sF 2 V TT 

2 


s 

- + — 


/•OO 

L„ (s' 


<5;.lo(^') 


Ml , iM, 

“NNLO ^ "NNLO ^2 


C^Wi _, Jh _ 

%NLOj^4 +%NL0^6 


Mk ' Jst, (s' - s - 4e)s''' 

s3 s^ /■<» ifiNyysx^ , 


■f 


j^NNLO 


(t) = 


b"° — + c"“ — 

''NNLOJ^f2_ +CNNLOJ^^4 

't" /•” ^O.LO(*') 

to {f - t- ie)tF 


l/2j_/,JV„ J_ , r„iv„ 

“ M|.l 




-NLOm4 


(s' - s - ie)s''‘ 

; r wj-°(f')^o%(t') ^^, 

Jto (f - t- ie)F 


r ^0°(*')^ofNLo(*') ,,, , *3 /°° jVo^"°(f')'?oyo(t') ^^, 


I Q “ “H — f 

^ Jto {f -t- ie)tF TT 

y°° ^o°(*')'^"lo(*') y^ ^ t'" j°°_ 

TT Jto {t'-t- ie)tF \ IT (t" -t' - ie)t"'^ 


{t' — t — 2e)i'" 

«i'Lo(<") 


1 dt' 


^NNLO(j( ^ , 


) + a 


iVi , ,1/2 


m3 


/Vf 


'w = 


^2 /-OO 

Ni ^ I 
"■NLO ^ m J. ■ \x,2 

TT (i' - i - ie)t' 

riVi ^ , 1/2 r, Ni _ 

"NNLO ^2 ^ "^1 ""NLO j^4 


-sy/oh') , t r Nr^it'WUt') 


■L 


it' — t — i€)t' 


dt\ 


SlUoit') 


TT {t' -t - ie)t'‘ 


,dt' 


Jto (i' 


' -syLoi*') 

-^- 2 '^ 

(s' - s - ie)s' 
ods' 




(t' — t — i€)t'‘ 


<2 /“ A“(<')-5y'/Lo(*')^^, 


(2 p ^NLO(t')5j_/^2^(i')^^^ ^ J2 ^00 ^^^(/'(^^//^(t') ^ i'2 


TT 


/ 


/■ 

•/io 


dNNLO 
it() 


w = - 


{t'-t- ie)t'^ TT {t'-t- ie)t'‘ 

3/2 t t^ 1“^ ^0^i'^')^o!ho(^') 


1 + 


)m/, 

{t" -t'- ie)t 


{f — t — ie)t‘ 
-^dt"]dt‘ 

7/2 
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f 

Jto 


" Ml 0 

' <3 /•“ <Lo(i') 


TT {f -t - ie)t'‘ 


73 

odt' I ( - 




to {t'-t - ie)t’' 


(3 p i?r°(t')5yyo(‘')^^, _ <! r ■ro°(<')Co(*') /, ^ r 


L 


^NNLO 

^NNLO 


(t' — t — ze)t'" 


Jto (f 


{f — t — ie)t'^ 


V ^ Jto it" 


{f — t — ie)t''^ 

it") 


„3/2 


- ie)t" 


odt" ]dt' 


(E.46) 


(E.47) 


where and SbJJl^Q are given by (147) and the remaining subtraction constants denote the quantities after 

the gauge transformation. Note that to and Co appear only in polynomial terms. In Mq, Mi and Mi, they can 
be reabsorbed into the NNLO subtraction constants. However, this is not the case for Nq, Ni, Ni and Rq. 
Here, we are required to fix the w-terms by imposing that the chirally expanded Omnes representation agrees 
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with the standard dispersive representation (or finally the two-loop representation). This somewhat awkward 
situation is just another manifestation of the fact that we identify the chiral representation with the Omnes 
dispersion relation although the phase shifts of the former have a wrong asymptotic behaviour. 

The comparison of the Taylor expansions of (E.47) and (69) leads to the relation (149) for the subtraction 
constants. 
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